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1 Introduction 



1.1 Background and motivation 



Problems in combinatorial enumeration and discrete probability can often be attacked by means of 
generating functions. If one is lucky enough to obtain a cloosed form generating function, then the 
asymptotic enumeration formula, or probabilistic limit theorem is often not far behind. Recently, 
several problems have arisen to which can be associated very nice generating functions, in fact ra- 
tional functions of several variables, but for which asymptotic estimates have not followed (although 
formulae were found in some cases by other means). These problems include random tilings (the 
so-called Aztec and Diabolo tilings) and other statistical mechanical ensembles (cube groves) as well 
as some enumerative and graph theoretic problems discussed later in the paper. 



A series of recent papers 



PW02 



PW04i IBP04| | provides a method for asymptotic evaluation of 



the coefficients of multivariate generating functions. To describe the scope of this previous work, we 
set up some notation that will be in force for the rest of this article. Throughout, we will assume 
that the generating function converges in a domain defining there a quasirational function 

P(Z) 



F(Z) 



Q(z)^n -=1^.(2)"^ 



a^Z'' 



(1.1) 



with polynomial P, Q, affine-linear i?j's, integer rij's and real s. (Here the quatities in boldface are 
vectors of dimension d and the notation Z'' is used to denote 0^=1 ^j^-) dimension three and 
below, we use X,Y,Z to denote ^1,^2 and respectively. We let V := {Z : g(Z) = 0} C 
denote the pole variety of F, that is, the complex algebraic variety where Q vanishes (Q will always 
be a polynomial) . Analytic methods for recovering asymptotics of ar from F always begin with the 
multivariate Cauchy integral formula 



{2mY 



Z-'"F- 



dZ 



(1.2) 



Here T is a d-torus, i.e. the product of circles about the origin in each coordinate axis (importantly, 
choice of a torus affects the corresponding Laurent expansion (|l.ip ). The pole set V is of central 
importance because the contour T may be deformed without affecting the integral as long as one 
avoids places where the integrand is singular. 

When V is smooth or has singularities of self-intersection type (where locally V is the union 
of smooth divisors), a substantial amount is known. The case where V is smooth is analyzed 



PW02]; th e exist ence under further hypotheses of a local central limit theorem dates back at 



least as far as 



BR83| . The more general case where the singular points of V are all unions of smooth 



components with normal intersections is analyzed using explicit c hanges of variables in PW04 | 
and by multivariate residues in BPOJ], the pre-cursor to which is Applicati ons in which 

V satisfies these conditions are abundant, and a number of examples are worked in [PW07|. For 
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bivariate generating functions, all rational generating functions fall we have seen fall within this class. 
Other local geometries are possible, namely those of irreducible monomial curve, e.g., — = 0, 
but, as will become clear, they cannot contribute to the asymptotic expansions, being non-hyperbolic. 

In dimension three and above, there are many further possibilities. The simplest case not handled 
by previous techniques is that of isolated quadratic singularity. The purpose of this paper is to 
address this type of generating function. All the examples in the title and in Section [4] are of this 
type. In fact, all the rational 3-variable generating functions we know of, that are not one of the 
types previous analyzed, have isolated, usually quadratic, singularities. The simplest case is when 
the denominator is irreducible and its variety has a single, isolated quadratic singularity; a concrete 
example in dimension d = 3 is the cube grove creation generating function, whose denominator 
= 1 + XYZ - {X + Y + Z + XY + YZ + ZX)/3 has the zero set illustrated in figure [TJ 



The main results of this paper. Theorems 13.71 and 13.91 below, are asymptotic formulae for the 
coefficients of a generating function having a divisor with this geometry. In addition to one or more 
isolated quadratic singularities, our most general results allow Q to be taken to an arbitary real 
power and we allow the possibility of other smooth divisors passing through the singluarities of Q. 
These generalizations complicate the exposition somewhat but are necessary to handle some of the 
motivating examples. 

As a preview of the behavior of the coefficients, consider the case d — 3 and F = 1/Q illustrated 
in figure [TJ The leading homogeneous term of Q in the variables {x,y,z) = (log X, log F, log Z) is 
xy + xz + yz. The outward normal cone (dual cone) at this point is the cone N* on which Q* > 0, 
where Q*{r,s,t) = {r + s + t)'^ — 2{r'^ + s"^ + t'^) is the dual quadratic to Q (see Section 12.61 for 
definitions). The asymptotics for Or in this example are given by Corollarv l3.8l for r in the interior 




Figure 1: an isolated quadratic singularity 
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of N* and by an easier result (Proposition [^221) when r ^ N*: 

f CQ*(r)-i/2 if J. is the interior of N* 

I exponentially small if r ^ N* 

The behavior of near 9N* is more complicated and is not dealt with in this paper. The generating 
function l/Q is the creation rate generating function for cube groves, discussed in Section [42l The 
edge placement generating function for cube groves (edge placement probabilities have more direct 
interpretations than do creation rates) has an extra factor of {1 — Z) in the denominator. Theorem l3.9l 
gives the asymptotics in this case, for r interior to N*, as 

Gr ^ C arctan6'(r) 

where is a homogeneous degree function of r whch can be expressed in terms of dual quadratic 
form Q* . Homogeneity of 9 implies that that there is a limit theorem a\r ^ 9{r) as A ^ cxd, where 
f :— r/|r| is the unit vector in the direction r. 

A total of five motivating applications will be discussed in detail in Section [H All of these may 
be seen to have factors with isolated quadratic singularities. There are known trivariate rational 
generating funtions with isolated singularities that are not quadratic. For example, the diabolo or 
fortress tiling ensemble has an isolated quartic singularity. Some of our results apply to this case, 
but a detailed analysis will be left for another paper. The last example goes slightly beyond what 
we do in this paper, but we include it because the analysis follows largely the same methods. 



Aztec diamond placement probability generating function |jPS98t 



F{X,Y,Z) 



Z/2 



(1 - YZ){i - f (X + AT-i + r + y-i) + Z2) 



(1.3) 



Cube groves edge probability generating function IpsosI ] 

2Z^ 



F{X, Y, Z) 



3(1 - Z){1 + XYZ - UX + Y + Z + XY + XZ + Y Z)) 



(1.4) 



Quantum random walk space-time probability generating function (aBN+01 ; BBBP08 | 

XZ - {1 + XY)Z'^ + YZ^) 



FiX,Y,Z) = 



(1 - Z2)(i ~{x + a:-i + Y + y-i)^/2 + Z2) 



(1.5) 
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Priedrichs-Lewy-Szego graph polynomial 



F{X, Y, Z) = [(1 - X){1 -Y) + {l-X){l-Z) + {l- Y){1 - Z)]-P (1.6) 



Super-ballot number generating function ( Ges92j | 



1 — IX 1 



1.2 Methods and organization 

Our methods of analysis owe a great debt to two bod ies of ex isting theory. Our approach to 



harmonic analysis of cones is fashioned after the work of jABG7C| . We not only quote t heir re sults 



on generalized Fourier transforms, which da te back somewhat farther to computations of Rie49l | and 



generalized function theory as decribed in |gS64| - but we also employ their results on hyperbolic 
polynomials to produce homotopies of various contours. Secondly our understa nding of the existence 



of these homotopies has been greatly informed by Morse theoretic results of |GM88|. We do not 
quote these results directly because our setting does not satisfy all their hypotheses, but the idea to 
piece together deforma tions local to strata is really the central idea behind stratified Morse theory 



as explained in [GM88j : see also the discussion of stratified critical points in Section [275 



An outline our methods is as follows. The chain of integration in the multivariate Cauchy 
integral (|1.2[) is a d-dimcnsional torus, T embedded in the complex torus C[J, where C, := C — {0}. 
Changing variables by Zj = exp(zj), the chain of integration becomes a chain C, the set of points 
with a fixed real part. Under this change of variables, the Cauchy integral ()1.2|) becomes 



^exp(-r.z)/(z)dz (1.8) 

where f := F o exp. Morse theoretic considerations tell us we can deform the chain of integration 
so that it is supported by the region where e^'^'^ is small (for large |r|) except near certain critical 
points. To elaborate, we can accomplish most of the deformation by moving x. The allowable region 
for such deformations of x is a component of the complement to amoeba of F (see Section 12.11 for 
definitions) . Heuristically we move x to the support point Xmax on the boundary of this region for a 
hyperplane orthogonal to r (see Sections 12.51 and those preceding for details) . Unfortunately, when 
Xmax is on the boundary (|1.8p fails to be integrable. Ignoring this, however, and continuing with 
the heuristic, we let q := Q o exp and hj := Hj o exp and we denote the leading homogeneous parts 
of q and hj by q and hj respectively. We then express / near x^ax as a series in negative powers of 
q and hj (this is carried out in Section 12. 7[) . Integrating term by term, each integral has the form 



exp(— r • z) dz 

c ' g(z)-/i(z)" 
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where /i" :— 11^=1 ^j^- Replacing z by iz, we recognize the Fourier transform of a product of a 
monomial with inverse powers of quadratics and linear functions. 

Much of what has been described thus far is based on known methods and results, most of which 
are collected in the preliminary Section [21 The bulk of the work, however, is in making rigorous 
these identities which involve Fourier integrals that do not converge, taken over regions which are 
not obviously deformations of each other (the part above where we said, "ignoring this, ..."). For 
this purpose, some carefully chosen d eformations are constructed, based largely on deformations 
found in Sections 5 and 6 of ABG70|. S pecificall y, we use results on hyperbolic polynomials (see 



Section 12.31 for definitions) established in jABG7Cll ] and els ewhere, to construct certain vector fields 
on C^. These vector fields, based on the construction of ABG70l . Section 5] and described in our 
Section [5.1[ then allow us to construct deformations in Section [5.2[ which satisfy several properties. 
First, they enact what Morse theory has guaranteed: they push the chain of integration to where 
the integrand of (|1.8p is very small, except near critical points, as in figure [5] below. Secondly, they 
do this without intersecting V, thereby allowing the integral to remain the same. This localizes the 
integral to the critical points, and allows us to concentrate on one critical point at a time. The 
resulting chain of integration is depicted in figure [21 





Real part = x 

^ max 

(a) in log space (b) in the original domain 

Figure 2: the localizing chain, in logarithmic and original coordinates 



Thirdly, they allow us to "straighten out" the chain of integration. Figure [31 shows that the 
chains in figure [21 as well as the original chain, are homotopic near the critical point to a (slightly 
perturbed) conical chain. 

Combined with the series expansion by homogeneous functions, this reduces all necessary inte- 
grals to a small class of Fourier -type integrals. Many of these are evaluated as generalized func- 

I and I 



tions m 



ABG70 



Rie49 



GS64| and elsewhere. In Section 



we summarize the relevant 



facts about generalized functions. The above deformations allow us to show, in Sections [ 
that these generalized functions, defined as integrals over the straight contour on the left of figure[3l 
do approximate the integrals we are interested in, which we must evaluate over the chains shown on 
the right of figure [Hand on figure [31 in order for the localizations to remain valid. Not all of the com- 
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Real part = u 



Figure 3: the projective chain 



putations we need are available in the literature. In Section TG-Gl we use a construction from |ABG70|, 
the Leray cycle, along with a resiude computation, to reduce the Foruier transform of \/{q ■ h) 
to an explicitly computable one-dimensional integral. It is this computation that is responsible for 
the explicit asymptotic formula for placement probabilities in the Aztec Diamond and Cube Grove 
problems. 

To summarize, the organization of the rest of the paper is as follows. Section [5] defines some 
notation in use throughout the paper, and collects preliminary results on amoebas, convex duals, 
hyperbolicity, and expansions by powers of homogeneous polynomials. Section [3] states the main 
results. Section [4] has five subsections, each discussing one of the five examples. The next two 
sections are concerned with the proofs of the main results. Section [5] constructs homotopies that 
shift contours of integration, while Section [6] evaluates several classes of integrals via the theory of 
generalized Fourier transforms. Finally, Section [7] concludes with a discussion of open problems and 
further research directions. 



2 Notation and preliminaries 



Several notions arising repeatedly in this paper are the logarithmic change of variables, duality 
between r and z, and the leading homogeneous part of a function. We employ some meta-notation 
designed for ease of keeping track of these. We use upper case letters for variables and functions in 
the complex torus C^, and lower case letters in the logarithmic coordinates. We will never use the 
notations without defining them, but knowing that, for example, F o exp will always be denoted / 
and z will always be log Z, should help the reader quickly recognize the setting. We will always use 
X and y for the real and imaginary parts of the vector z g C*. Boldface is reserved for vectors. The 
leading homogeneous part of a function is denoted with a bar. Rather than considering the index r 
of flr to be an element of Z , we consider it to be an element of a space (R'')* that is dual to the 
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domain in which z hves, with respect to the pairing r • z (the space (R'')* is a subset of the full 
dual space (C*)* but all our dual vectors will be real). Many functions of r use in what follows are 
homogeneous degree 0; letting f denote the unit vector r/|r| we will often write these as functions 
of r. The logarithm and exponential functions are extended to act coordinatewise on vectors. Thus 

exp(zi, . . . , Zd) := (cxp(zi), . . . ,exp(zrf)) ; 
log(Zi,...,Z<i) := (log(Zi),...,log(Zrf)). 

We also employ the slightly clunky notation 

ReLogZ (log |Zi|, . . . , log = Re{logZ} 



for the coordinatewise log-modulus map, having found that the notations in use in GKZ94| do not 
allow for quick visual distinction between log and ReLog . 

Our chief concern is with generating functions that are ordinary power series, convergent on the 
unit polydisk, and whose denominator is the product of smooth and quadratically singular factors 
that intersect the closed but not the open unit polydisk. It costs little, however, and there is some 
benefit to work in the greater generality of Laurent series representing functions with polynomial 
denominators. Indeed, Laurent series arise naturally in the examples (though these Laurent series 
have exponents in proper cones, and may therefore be reduced by log-affine changes of coordinates 
to Taylor series). 

Definition 2.1 (homogeneous part). For analytic germ f : (C'^,z) — > C at a point z G C^, we let 
deg(/, z) denote the degree of vanishing of f at z. This is zero if f(z,) ^ and in general is the 
greatest integer n such that /(z + w) — 0(|w|") as w ^ 0. Also, deg(/, z) is the least degree of any 
monomial in the ordinary power series expansion of f{z + •) around 0. We let flhom(f,2,) denote 
the sum of all monomials of minimal degree in the power series for /(z + •) and we call this the 
homogeneous part of f at z. Thus 

/(z + w) = flhomif, z)(w) + O (\wf°s(f,^)+i^ 
for small \w\.When z — 0, we may omit z from the notation: thus, phom{f) :— flhom{f,Q). 

A number of connections between zeros of a Laurent polynomial F, the Laurent series for l/F, 
the Newton polytope for F and certain dual cones to this polygon were worke d out in t he 1990's by 



Gelfand, Kapranov and Zelivinsky. We summarize some relevant results from GKZ94 Chapter 6] 



2.1 The Log map and amoebas 



Let F be a Laurent polynomial in d variables. Let Vf denote its zero set in C^. Let (C*)'' C 
denote the d-tuples of nonvanishing complex numbers. Following GKZ94| we define the amoeba 
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of F to be the image under ReLog of the zero set of F: 

fianioeba(F) {ReLog z : z e Vf n (C*)'*} C 



The simplest example is the amoeba of a linear function, such as F — 2 — X — Y, shown in figure 4(a) 
The amoeba of a product is the union of amoebas, as shown in figure [4(b) [ 




3 
2 






4 -3 -2 -i U 

-2 
-3 


utlJI / 2 3 4 



(a) Bamocba(2 - X - Y) 



(b) Bamocba(3 - 2X - Y)(3 
X - 2Y) 



Figure 4: two amoebae 



The rational function 1/F has, in general, a number of Laurent se ries expa nsions, each convergent 
on a different subset of C*. Combining Corollary 1.6 in Chapter 6 of jGKZ94| with Cauchy's integral 
theorem, we have the following result. 

Proposition 2.2. The connected components ofW^\fiamoeha{F) are convex open sets. The com- 
ponents are in bijective correspondence with Laurent series expansions for 1/F as follows. Given 
a Laurent series expansion of 1/F, its open domain of convergence is precisely ReLog where 
B is a component of \ Pamoeha[F) . Conversely, given such a component B, a Laurent series 
1/F = ^flrZ'' convergent on B may be computed by the formula 



1 



1 



■dZ 



{2mY Jt~ F{Z) 

where T is the torus ReLog ~^(x) for any x G _B. Changing variables to Z ~ exp(z) and dZ = Zdz 
gives 

^ ^ — ^ dz (2.1) 

□ 



1 



where f = F o exp and Tr is the torus i?'*/(27rZ)'*. 



2.2 Dual cones, tangent cones and normal cones 

Let (R'^y denote the dual space to R'* and for y e M'' and r e (R'')*, use the notation r -y to denote 
the pairing. Let L be any convex open cone in R'^. The (closed) convex dual cone L* C (R'^)* is 
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defined to be the set of vectors v e (M'')* such that v • x > for all x e L. Familiar properties of 
the dual cone are 

L C M L* DM*; (2.2) 

{LnM)* = ch(L*UM*). (2.3) 

Suppose X is a point on the boundary of a convex set C. Then the intersection of all halfspaces that 
contain C and have x on their boundary is a closed convex afRne cone with vertex x (a translation 

by X of a closed convex cone in M."^) that contains C. Translating by — x and taking the interior 
gives the (open) tangent cone to C at x, denoted by i3tanx(C). An alternative definition is 

6tanx(C) = {v : X + ev e C for all sufficiently small e > 0} 

(where B is the unit ball). The (closed) normal cone to C at x, denoted N*x(C), is the convex 
dual cone to the negative of the tangent cone: 

N*x(C) = (-6tanx(C))* . 

Equivalently, it corresponds to the set of linear functionals on C that are maximized at x, or to the 
set of outward normals to support hyperplanes to C at x. 

Definition 2.3 (proper dual direction). Given a convex set C, say that r is a proper direction for 
C if the maximum ofr-x on L is achieved at a unique Xmax S C. We call Xmax the dual point for 
r. The set of directions r for which r - x is bounded on C hut r is not proper has measure zero. 

The term tangent cone has a different meaning in algebraic contexts, which wc will require these 
as well. (The term normal cone has an algebraic meaning as well, which wc will not need.) To avoid 
confusion, we define the algebraic tangent cone of / at x to be V^Yiomu z) ■ equivalent but 
more geometric definition is that the algebraic tangent cone is the union of lines through x that are 
the limits of secant linos through x; thus for a unit vector u, the line x+iu is in the algebraic tangent 
cone if there arc x„ G V/ distinct from but converging to x for which (x„ — x)/||x„ — x|| ^ ±u. 
This equivalence and more is contained in the following results. We let denote the unit sphere 
{(zi, . . . ,Zd) ■ l^ip H + = 1} and let Sr ■= rSi denote the sphere of radius r. 

Lemma 2.4 (algebraic tangent cone is the limiting secant cone). Let q be a polynomial vanishing 
to degree m>l at the origin and let q := fihom,{q) be its homogeneous part; in particular, 

q{z) = q{z) + R{z) 

where q is a nonzero homogeneous polynomial of degree m and R{z) = 0(|z|'"+^). Let denote the 
polynomial 

q,{z) := e-'"(7(ez) =q(z) + i?,(z) 

where Re{z) = e~'"i?(ez) ^0 as e ^ 0. Let Ve '■= Vq^ fl 5*1 denote the intersection of {q^ = 0} with 
the unit sphere. Then Ve converges in the Hausdorff metric as e — > toVqCi Si. 



9 



Proof: On any compact set, in particular 5*1, ^ uniformly. If z*^"^ z and z*^"^ G Vi/„ then 
for each n, 

|g(z("))| = |gi/n(z("^) + i?i/„(z("))| - |i?i/„(z("))| ^0. 

Hence g(z) = by continuity of q and we see that any limit point of Ve as e is in Vq H Si. 
Conversely, fix a unit vector z G Vq. The homogenous polynomial q is not identically zero, therefore 
there is a projective line through z along which q has a zero of finite order at z. Back in affine space, 
there is a complex curve 7 in the unit sphere along which q is holomorphic with a zero of some finite 
order fc at z. As e ^ 0, the holomorphic function goes to zero uniformly in a neighborhood of z 
in 7; hence there are k zeros of converging to z as e — > 0, and therefore z is a limit point of as 
e ^ 0. □ 



2.3 Hyperbolicity for homogeneous polynomials 



The notion of hyperbolic polynomials arose first in GarSd ] in connection with solutions to wave-like 



partial differential equati ons. The same property turns out to be very important as well for convex 



programming, cf. Giil97j from which much of the next several paragraphs is drawn. 



Let / be a complex polynomial in d-variables and f{D) denote the corresponding linear partial 
differential operator with constant coefficients, obtained by replacing each Xj by d/dxj . For example, 
if / is the standard symplectic quadratic iS'(x) = x\ — '}2'^i=2 ^? ^'^^'^ ^i^) is the wave operator 
{d/dxiY — J2j=2i^/^^j)^- Garding's object was to determine when the equation 

f{D)u = g (2.4) 

with g supported on a halfspace has a solution supported in the same halfspace. The wave operator 
has this property, and in fact there is a unique such solution for any such g. It turns out that the class 
of / such that ()2.4p always has a solution supported on the halfspace is precisely character ized by 



Hor9C 



the property of being hyperbolic, as defined by Garding. In this case, it was later shown ( 
Theorem 12.4.2] that the solution is in fact unique. The theory of hyperbolic polynomials serves in 
the present paper to prove the existence of deformations of chains of integration past points of the 
pole manifold at which the pole polynomial is locally hyperbolic, in order to 

We begin with hyperbolicity for homogeneous polynomials, which is a simpler and better devel- 
oped theory. As before, we use A rather than / for a homogeneous polynomial. 

Definition 2.5 (hyperbolicity). Say that a homogeneous complex polynomial A of degree m > 1 is 
hyperbolic in direction v G M.'^ if A{v) ^ and the polynomial A{x + tv) has only real roots when x 
is real. In other words, every line in real space parallel to v intersects Va exactly m times (counting 
multiplicities). 

While seemingly weaker, the requirement of avoiding purely imaginary roots is in fact easily seen 
to be equivalent. 
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Proposition 2.6. Hyperbolicity of the homogeneous polynomial A in the direction v is equivalent 
to the condition that v4(v + iy) ^ for all y G R'*. 

Proof: Because A is homogeneous, when A 7^ 0, we have A{Xz) if and only if A{z) — 0. With 
A = i • s, a purely imaginary number not equal to zero, we see that A{v + iy) 7^ for all y G M.'^ is 
equivalent to A{y + isv) 7^ for all y G K'' and all nonzero real s. Writing z — t + is, this becomes 
^(y + tu + isu) for all y G M'' and real s 7^ 0, which is equivalent to A{y + zu) ^ when z is 
not real, which is the definition of hyperbolicity. □ 



The further properties we need are well known and are proved, among other places, in Giil97l . 
Theorem 3.1]. 

Proposition 2.7. The set of v for which A is hyperbolic in direction v is an open set whose 
compenents are convex cones. Denote by 'K.^(A) the connected component of this cone that contains 
a given v. Some multiple of A is positive on 'K'^{A) and vanishing on dT^^{A), and for x. G K'^(yl), 
the roots of A(x + tv) will all be negative. □ 

Semi-continuity properties for cones of hyperbolicity play a large role in the construction of 
deformations. A lower semi-continuous function / satisfies f{x) < liminfy^j; f{y). The property is 
important in elementary analysis because a lower semi-continuous function on a compact set achieves 
its infimum; generalizing to set-valued functions, the conclusion is roughly that the empty set is not 
a limit value and therefore that a continuous section can be defined. In this section, we develo p 



semi-continuity properties for cones of hyperbolicity (a topic that occupies many pages of ABG70|). 



The following proposition and definition define a family of cones {K"^'^(x)}xgR<j which will be 
used to prove two critical semi-continuity results for cones of hyperbolicity for log-Laurent polyno- 
mials (Theorem 12 . 141 below) . 

Proposition 2.8 (first semi-continuity result). Let A be any hyperbolic homogeneous polynomial, 
and let m be its degree. Fix x with A(x) = and let a := flhom{A, x) denote the leading homogeneous 
part of A at x. If A is hyperbolic in direction u then a is also hyperbolic in direction u. Consequently, 
if B is any cone of hyperbolicity for A then there is some cone of hyperbolicity for a containing B. 



Proof: This follows from the conclusion (3.45) of |ABG7d . Lemma 3.42]. Because the development 
there is long and complicated, we give here a short, self-contained proof, provided by J. Borcea 
(personal communication). If P is a polynomial whose degree at zero is k, we may recover its 
leading homogeneous part fihom(P) by 

i3hom(P)(y) = lim A''P(A- V) • 

A — >oo 

The limit is uniform as y varies over compact sets. Indeed, monomials of degree k are invariant 
under the scaling on the right-hand side, while monomials of degree k -\- j scale be A~^ , uniformly 
over compact sets. 
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Apply this with P(-) = A{x + •) and y + in place of y to see that for fixed x, y and u, 

a(y + tu)= lim A''A(x + A^^ (y + tu)) 

A — ^oo 

uniformly as t varies over compact sub-intervals of R. Because A is hyperbolic in direction u, for 
any fi xed A, all the zeros of this polynomial in t are real. Hurwitz' theorem on the continuity of 
zeros Con78l , Corollary 2.6] says that a limit, uniformly on bounded intervals, of polynomials having 
all real zeros will either have all real zeros or vanish identically. The limit, a(y + tu) has degree 
k > 1; it does not vanish identically and therefore it has all real zeros. This shows a to be hyperbolic 
in direction u. □ 

Definition 2.9 (family of cones in the homogeneous case). Let A be a hyperbolic homogeneous 
polynomial, let B he a cone of hyperbolicity for A. If A{x) — 0, define 

K^'^(x) 

to be the cone of hyperbolicity of flhom{A,x.) containing B, whose existence we have just proved. If 
A{x) y^Owe define K^'^{x) to be all ofR'^. 

As an example of a hyperbolic polynomial, let S ~ J be the standard Lorentian 

quadratic. Then K^^^ {S) is the Lorentz cone {x : xi > -^/xl + • • • + a;^}. Any quadratic of Lorentizan 
signature is obtained from this one by a real linear transformation; we see therefore that for any 
Lorentzian quadratic, the boundary of the cone of hyperbolicity is the algebraic tangent cone. 

The class of hyperbolic polynomials in a given direction direction v is closed under products, 
and K"(AA') — K"(A) n K"(A'). The class contains all linear polynomials not annihilating v and 
all real quadratic polynomials p of Lorentzian signature for which p(v) > (v is time- like) . 



2.4 Hyperbolicity and semi-continuity for log-Laurent polynomials on the 
amoeba boundary 

For a function that is not locally homogeneous, there are two natural generalizations of the definition 
of hyperbolicity. Both are equivalent to the notion of hyperbolicity already defined, in the case of a 
homogeneous polynomial. Useful features of the two definitions are revealed in the subsequent two 
propositions. 

Definition 2.10. Let / ; — > C vanish at z and he holomorphic in a neighborhood of z. We 
say that f is strongly hyperbolic at z in direction of the unit vector v if there is an e > such 
that /(z + tv' + iu) 7^ for all real < t < e, all |v' — v| < e, and all u G M'' of magnitude at 
most e. In this case we may say that f is strongly hyperbolic at z in direction v with radius e. 
Say that f is weakly hyperbolic in direction v if for every M > 0, there is an e > such that 
/(z + tv + iu) for all real < t|v| < e, and for all u G R'' of magnitude at most e additionally 
satisfying |u|/(t|v|) < M. 
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Proposition 2.11. Let A = fihom{f,z). Then A is hyperbolic in direction u if and only if f is 
weakly hyperbolic in direction u at z. 

Proof: The homogeneous polynomial A fails to be hyperbolic at in direction u if and only if there 
is some real y such that ^(u + iy) = 0. By Lemma this happens if and only if /(z + w„) = for 
some sequence {w„} converging to with w„/w„ converging to (u + iy)/|u + iy|. This is equivalent 
to failure of weak hyperbolicity of / at z in direction u. □ 

Remark. It is immediate from the definition that strong hyperbolicity is a neighborhood property: 
if / is strongly hyperbolic at x + zy in direction v with radius e, then for |y' — y| < e and |v' — v| < e, 
/ is strongly hyperbolic at x + iy' in direction v' with direction e — max{|y' — y|, |v' — v|}. Weak 
hyperbolicity of / at z in direction v extends to a neighborhood of v by Propositions 12.71 and 12.111 
Extending weak hyperbolicity to neighboring z is much trickier. 

Proposition 2.12. Let F be a Laurent polynomial in d-variables. Suppose that B is a component of 
flamoeba(F) and x G dB. so that f := Focxp vanishes at some point x+zy. Let f := flhom{f, x+zy) 
denote the leading homogeneous part of /(x + iy + ■)■ Then f is strongly hyperbolic at x + iy, some 
complex scalar multiple of f is real and hyperbolic, and some cone of hyperbolicity K"(/) contains 
ptan^iB). 

Proof: Strong hyperbolicity of / in any direction u e fitanx(i?) follows from the definition of 
the amoeba. Strong hyperbolicity is stronger than weak hyperbolicity, hence hyperbolicity of / in 
direction u follows from Proposition 12.111 The vector u G i3tanx(i?) is arbitrary, whence K"(/) 3 
6tanx(i?). To see that some multiple of / is real, let u be any real vector in fitanx(i3), let m denote 
the degree of /, and let 7 denote the coefficient of the term of A{fu + y). Then 7 is the degree 
TO coefficient of f{zu), hence is nonzero and does not depend on y. For any fixed y, the fact that 
/(zu + y) has all real roots implies that the monic polynomial 7^^/(zu + y) has all real coefficients. 
□ 

Definition 2.13 (hyperbolicity and normal cones at a point of V/). Let F be a Laurent polynomial, 
B a component of \ fiamoeba{F), and Z = exp(x + iy) G V/ with x G dB. We let f :^ F o exp 
and let 

K/'^(Z) :=K"(i?/iom(/,x + zy)), (2.5) 

denote the (open) cone of hyperbolicity of f := fihom{ f,:x. + iy) that contains B, whose existence is 
guaranteed by Provosition \2. ISl We also define the normal cone 

N*(Z) := (N*)^'-^(Z) = (Kf^^iZ))* . (2.6) 

We see immediately from Proposition 12.121 that 

K^^-^(Z) D tanx(S) (2.7) 
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and hence 



N*(Z) C N\{B). 



In order to produce deformations, we will need to know that the cones vary semi- 

continuously as Z varies over the torus exp(x + iW^). We have seen already that all of these cones 
contain 13tanx(i3). What is needed, therefore, is an argument showing that K-^'^(Z') contains any 
u g K-^'^(Z) when Z' is sufficiently close to Z and u ^ fitanx(i3). In fact, not every polynomial 
admits a semi-continuous choice of tangent subcone; a counterexample is xy + z^. However, in the 
case where x G dB, we are able to use strong hyperbolicity in directions v e fitanx(-B) to prove 
semi-continuity even outside of f5tanx(-B). The main result of this section is exactly such an analogue 
of Proposition 

Theorem 2.14. Suppose that an analytic function f is strongly hyperbolic in direction v at the 
point z = x + iy. Let f := phom{f ,z) . Let u £ K^(/) be any point in the cone of hyperbolicity of 
f containing v. Then f is strongly hyperbolic in direction av + (1 — a)u for every < a < 1. 

Corollary 2.15. 



(i) If B is a component of flamoeba{FY and x e dB, then K.f'^{Zi) is semi- continuous in Z as 
y varies with Z = exp(x + iy), meaning that K-^'^(Z) C liminfz' K.-^'^{Z'). 

{a) If A is a homogeneous polynomial and B is a cone of hyperbolicity for A, then K"^'^(y) is 
semi- continuous in y. 



Proof: Pick any v e 13tanx(i?). The function / is strongly hyperbolic in direction v, hence by 
Theorem 12.141 it is strongly hyperbolic at x + iy in every direction u e K'^(/). Because strong 
hyperbolicity is a neighborhood property, it follows that for every y' in some neighborhood of y, 
some cone of hyperbolicity of / contains K'^(/). All these cones contain v, hence these are the 
cones K-'''^(Z') (with Z' := exp(x -I- iy')), and hence all these cones contain K-^'^(Z). The proof in 
the homogeneous case is analogous, again because each a fihom(^,y) is strongly hyperbolic in 
direction v for any v G B. □ 



Proof of Theorem 12. 141 The pr oof is based on a technique of Carding [GarSOl . Theorem H 5.4.4] 
that is used in the proof of ABG70l . Lemma 3.22]. Let / be strongly hyperbolic at x + iy in direction 
V with radius e and choose any u G K'^(/). For the remainder of this argument, we assume that y' 
and v' satisfy 



V - V < 



a consequence is that / is strongly hyperbolic in direction v' at x -I- iy' with radius e/2. For any 
b G M'^, if s is purely imaginary with |s||b| < e/2, then the imaginary vector sb -I- z(y' — y) will have 
magnitude less than e. By hypothesis, when < < < e, the function 



si^ f{y. + iy' + sh + t\') 



(2.8) 



14 



will therefore be nonzero. 

As the complex argument s tends to zero, there is an expansion 

/(x + iy + s(av' + (1 - a)u)) = s"7(av' + (1 - a)u) + s"'+^B{a, s) 

where B is analytic. The homogeneous function / does not vanish on the convex hull of u and the 
(e/2)-ball about v, hence \f{av' + (1 — a)u)\ is uniformly bounded away from zero for a € [0, 1] and 
|v' — v| < e/2. It follows that for a sufficiently small S (which we take also to be less than e), the 
function 

s t-^ /(x + iy' + s{av + (1 - a)u) + tv) 

has exactly m roots bounded in absolute value by 5, as long as |y' — y| and t are both boimded in 
magnitude by S. Once 2(5|av' + (l — a)u| < e for all < a < 1, then, taking b — q:v' + (1— q:)u in (|2.8p . 
we see that these m roots cannot be purely imaginary, and their real parts must therefore retain the 
same sign as a, (3 and y' vary. When a = 1, these are the m roots in s of /(x + iy' + (s + i)v'), so 
the real parts are t less than the real parts of the roots of /(x + iy' + sv') which are all negative by 
strong hyperbolicity of / at x + iy' in direction v'. We conclude that for all positive real s in the 
interval < s < (5, the function /(x + iy' + s(q;v' + (1 — a)u)) does not vanish, finishing the proof 
of strong hyperbolicity with neighborhood size 5, for any a G [0, 1]. □ 

Corollary 2.16. Let F be a Laurent polynomial and f := i^oexp. Let x G dB for some component 
B of Pamoeha{FY . Let he a continuous unit section o/ K-'^^^(exp(x + i-)). In other words, 9 : 
(R/(27rZ))'' — > S'^^^ is continuous and 9{y) € K-^'"^(exp(x + iy)) for each y. Then there is some 
eo > such that for all < e < cq, /(x + iy + eO{y)) ^ 0. 

Proof: For each y, let e(y) be a radius of strong hyperbolicity for / at x + iy in direction 0(y). 
Choosing a neighborhood M(jf) such that \9{y') — 9(j)\ < e(y)/2 when y' e A/'(y), we see that 
e(y)/2 is a radius of strong hyperbolicity for / at x + iy' for any y' E A/'(y). Covering the 
compact set (]R/(27rZ))'^ with finitely many neig hborhoods 7V(y(i'), . . . ,7V(y(")), we may choose 
eo = min„e(y(")). □ 

Examples and counterexamples 

It is important to understand how semi-continuity may fall short of continuity. This is illustrated in 
the following examples. To avoid misleading you with the pictures, we note that all of the upcoming 
figures show complex lines in C^, but that for obvious dimensional reasons, only the intersection 
with the E X M subspace is shown. 

Example 2.17 (cones drop down on a substratum). Let F — L1L2 = {3 — X — 2Y)(3 — 2X~Y). This 
differs from figure\^in that now I2 also passes through (1, 1). However, since L2 in this example is the 
inversion of L2 in example \2.19[ the amoeba is the same as in example \2.19\ . We will see that the cone 
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Figure 5: the zero set of the function L1L2 := (3 - X - 2Y){i -2X-Y) 

K(Z) drops discontinuously as Z — )■ (1,1), in contrast to examvle \2.1iA The subset ofVp lying in 
ReLog^^B) is the union of two rays{{Y-l) = -{X-l)/2 : X < 0}U{(r-l) = -2{X-1) : X > 0} 
with the common endpoint (1,1). For any point Z in this set other than (1,1), the cone K(Z) is 
equal to taniogz^ which is a halfspace. For Z = (1, 1), the cone K(Z) is still equal to taniogz(^), 
but now this is a proper cone bounded by rays with slope —1/2 and —2. This cone is the intersection 
of the two halfspaces that are possible values of the cone at nearby points, thus K(l,l) is equal to 
the Urn inf o/K(Z) for nearby Z, but there is a discontinuity at (1, 1). 

Compare this to examvle \2.19\ Here, Vp H ReLog^^(_B) is the union of two rays with different 
endpoints (1,1) and (—1,-1) and K(Z) is continuous, being constant on each ray and equal to a 
different halfspace on each ray. 

The containment tanx(-B) C K'^'-^(Z) for Z ~ exp(x + iy) £ V/ may be strict. We will see later 
that this causes a headache, so we formulate a property allowing us to bypass this trouble in some 
cases. 

Definition 2.18. Say that x is a well covered point of dB if K^'^{Z) — tanx(-B) for some Z = 
exp(x + iy). 

We now give two examples of points that arc not well covered. 

Example 2.19 (two lines with ghost intersection). Let F = L1L2 = (3 - X - 2y)(3 + 2X + Y). 
The variety Vp is shown on the left of figure\d[ Its amoeba is identical to the amoeba on the right of 
figure^ Indeed, it is the union of flamoeba{S-' X ~2Y) and flamoeba{3 + 2X + Y) , the latter of which 
is identical to flamoeba{3 — 2X — Y) because the amoeba of F{—X,—Y) is the same as the amoeba 
of F{X,Y). The component B of M."^ \ fiamoeba{F) containing the negative quadrant corresponds to 
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L2 

(a) Vf 



(b) a component of Bamoeba(F) 



Figure 6: the zero set of (3 — X — 2y)(3 + 2X + Y) from example 12.191 and the OPS component 

the ordinary power series. An enlargement of this component is shown on the right of figure\^ For 
X 7^ (0, 0) G dB, the only point Z = exp(x + iy) of Vf is the real point Z = ± exp(x), the positive 
point being chosen for the part of dB in the second quadrant and the negative point for the part of 
dB in the fourth quadrant. In either case, K(Z) is equal to the half space tanx(i?). 

On the other hand, whence — (0,0), the linearization off at x. is just iii2 :— {X + 2Y)(2X + Y). 
The zero set of which contains the two rays forming the boundary of 

tanx(B) = {{u, u) e : 2u + w < and u + 2v <0} . 

There are two points Z G Vf in ReLog "^(0, 0), namely (1,1) and (—1,-1). The first is in Vli 
and the second is in Vl2- The cone K(l, 1) is the halfspace {(m, u) € : w + 2u < 0}, while the 
cone K(— 1, —1) is the halfspace {(u,w) € : 2w + t; < 0}. Both of these cones stricly contain the 
cone tanx(i?). The term "ghost intersection" refers to the fact that the two curves ReLog V^^ and 
ReLogVia intersect at (0,0) but the lines V^i and Vl2 have different imaginary parts and have no 
intersection on the unit torus (though they do intersect at (—3,3)^. 

Next we include an example which is the closest we can get in two dimensions to the amoeba of 
a quadratic point (which can occur only in dimensions three and higher). 

Example 2.20 (critical set has large intersection with a torus). Let F = 1 — ■\/l/2(l — X)Y — XY^ 
be the denominator for the generating function for a one- dimensional Hadamard quantum random 
walk (see (BPOl l). The component B ofM.'^\Pamoeba{F) corresponding to the ordinary power series 
is that component of the complement of the shaded region in figure which contains the negative 
quadrant. 

To illuminate this example a little more, observe that as we go around the boundary of the 
amoeba, starting at the origin and leaving to the northwest, the dual cone is a single projective 
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Figure 7: the amoeba for F = 1 - {1 - X)Y/V2 - XY'^ 



direction A € at every point other than the origin. Parametrizing KF'^ by X = y/x, we see A 
decreasing from c := (1— ^/l/2)/2 to zero as the tentacle goes to infinity, then from to — oo coming 
back down the other side of tentacle and from +oo to 1 going up and out the northwest tentacle, 
and so forth. For each point of d{fiamoeba{F)) other than the origin, there is a unique x G and 
y G with /(x + iy) = 0; the cone K(exp(x + iy)) is equal to the halfspace tanx(-B). 

On the other hand, when x = (0,0), the cone tanx(i?) is bounded by the two rays A = c and 
A = 1 — c. This is noted in figure by the arrow matching the interval [c, 1 — c] to the single 
point at the origin. It is easy to check that if {X,Y) G then \X\ — 1 if and only if \Y\ — 1. 
Thus the intersection of Vf with the unit torus is the smooth topological circle parametrized by 
{(,^(e^^),e^^) :2/GR/{27rZ)}. 

As {x,y) varies over this curve, the cone K(e*^,e*^) remains a halfspace, the slope of whose 
normal varies smoothly between (1 + •\/l/2)/2 and (1— ^l/2)/2 and back. All of these cones strictly 
contain tan^iB). Thus the cone tanx(i3) is the intersection of the cones {K(Z) : Z G Vf H T^} but 
these all strictly contain tanx(i?). 

2.5 Critical points 

It is time to give further examination to the role of x^ax- The modulus of the term Z""" in the Cauchy 
integral is constant over tori, and among all tori in ReLog ~^(-B) the infimum of |Z^''| occurs on the 
torus ReLog "^(xmax)- This already indicates that this torus is a good choice, but we may get some 
more intuition from Morse theory. The space V is a Whitney stratified space: a disjoint union 
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of smooth real manifolds, called strata, that fit together nicely. The axioms for this may be found 



in Section 1.2 of part I of GM88[, along with some consequences. We will use Morse theory only as 

(2.9) 



a guide, quoting precisely one well known result, namely local product structure: 

A point p in a /c-dimensional stratum of a stratified space V has a neighborhood 



in which V is homeomorphic to some product S x X . 
This is needed only for the proof of second part of Proposition 12.221 below, which i n turn is used 



only for classifiying critical points when computing examples. According to GM88j . a proof may 
be found in mimocgraphed notes of Mather from 1970; it is based on Thom's Isotopy Lemma which 
takes up fifty pages of the same mimeographed notes. 

A point Z e V is a critical point for the smooth function h ii dh\s vanishes at Z, where S is the 
stratum containing Z. Goresky and MacPherson show that in fact such points are the only possible 
topological obstacles to lowering the value of h. Taking h — — f • ReLogZ, we see that (z) if there 
is no critical point in ReLog "^(xmax) then this torus is in fact not the best chain of integration, 
and (a) if there is a critical point in this torus then we may use this fact to help us compute Xmax- 
Because we do not give a rigorous development of stratified Morse theory here, we give a definition 
of the critical set in terms of cones of hyperbolicity, then indicate the relation to Morse theory. 

Definition 2.21 (minimal critical points in direction r). Fix a Laurent polynomial F in d variables 
and a component B of o/M'* \fiamoeha{F). For a proper direction r, let Xijiax(f) denote the unique 
point on dB maximizing r • x and let Vi = Vi(f) = Vi(x,jiax) denote the intersection of V with 
ReLog "^(xniax)- Recall the notation N*(Z) for the dual cone to the cone K''''^(Z) and define the 
set of minimal critical points by 

ficrit{r) := {Z G Vi(r) : r € N*(Z)} . 

A logarithmic version of ficrit is 

flW{r) := {y e Tk : exp(x,„ax + iy) e flcrit{r)} . 



The te rm "minimal" refers to the fact that ReLog Z e dB and following the terminology of PW02 



PW04 | 



Proposition 2.22. Fix a Laurent polynomial F in d variables, let f :— F o exp, and let B be a 
component of \ fiamoeba{F) . If Z £ "^"1(1") is not in f!crit{r) then there is some v £ Ji.-^'^{Zi) with 
f • V = 1. Conversely, i/Z e ficrit{r) then Z is a critical point for the function (j) := r ■ logZ on the 
stratified space V. 



Proof: If Z ^ ficrit(r) then by definition of the dual cone, the maximum of f • x on fitanx(i?) is 
strictly positive. Letting v' denote a vector in i3tanx(i?) for which r • x > 0, we may take v to be 
the appropriate multiple of v'. 



19 



For the converse, suppose that Z is not a critical point of the function cj) on V. Then z := logZ 
is not a critical point for f :— F o exp on log V; denoting d{(j3 o exp) by r, we see, by definition of 
criticality in the stratified sense, that r\s is not identically zero, where S is the stratum of log V in 
which z lies. 



We claim that the linear space T'z(5') is what ABG7Cll | call a lineality for the function f° := 



6hom(/, z), meaning that /°(w + w') = /°(w) for any w' S T'z(S') and any w e C^. To see this, for 
any w G C^', let w = W[| +wj^ denote the decomposition into an element W|| G Tx(z) and an element 
in the complementary space Tx{z)^ . Write / as a power series ^ CrW|| in W|| with coefficients that 
are power series in wj^. The coefficients Cr(0) vanish for |r| < m := deg(/, z). By (|2.9p . the degree 
of vanishing of / at any point of S is the same, hence Cr(wx) vanish identically for r < m. This 
implies that the only degree m terms in the power series for / near z are those of degree m in W|| , 
which implies that /°(w) depends only on W||, proving the claim. 



By Proposition 12.121 we know that f° is hyperbolic. By |ABG70l . Lemma 3.52], the real part 
of the linear space Tx{S) is in the edge of cK^-'^{Z), meaning that K-^'^(Z) is invariant under 
translation by vectors in Tx{S). Any real hyperplane not containing the edge of a cone intersects 
the interior of the cone. Applying this to the real hyperplane {x : r • x = 0}, (recall by assumption 
of noncriticality that this hyperplane does not contain Tx{S)), we conclude that there is some point 
p G cK-f-^{Z) with r • p = 0. This implies Z ^ ficrit(r). □ 

Showing that 6crit(r) is contained in the set of critical points of the logarithmic gradient enables 
us to use algebraic computational methods, cf. the A ztec Diamond computations in Section 14.11 



Some of this algebraic apparatus is detailed further in jPW04l : |BP04| |: for the present purpose, the 
following observations will suffice. When Z is a smooth point of Vp, the homogeneous part / of 
f := F o exp is a linear map vanishing on the tangent space to / at logZ. Hence the cone of 
hyperbolicity of / is an open halfspace, and the dual is the normal vector to this halfspace, which 
is the logarithmic normal to Vf at Z. (Thus in some sense, the dual cone N* is a set-valued 
generalization of the logarthmic gradient map.) To compute the smooth points of ficrit(r), we 
observe that the gradient of r • logZ is (ri/Zi), . . . , {rd/Zd)- Thus, for Z to be a smooth critical 
point, on the divisor {Hj = 0} we must have 

H, = 0; 

Similarly, for a stratum which is the transverse intersection of k smooth divisors {Hj '■ 1 < j < k} 
with logarithmic normals ViogHj, the equations for critical points in direction r are Hi = ■ ■ ■ = 
Hk = and 

r G (Viogffi,...,Viogi/fc), (2.11) 

the linear span of the k logarithmic gradients. Generically, this defines a zero-dimensional variety, 
meaning that the number of solutions is finite and nonzero. 
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For functions / and 5, define the notation 



f = o,^^{g)^\f{x)\<e-^-g{x) 
for some /3 > and all sufficiently large x. 

Proposition 2.23. Let F be a Laurent polynomial and ^^.arZ"" he a Laurent series for 1/F, 
convergent on a domain ReLog^^(i?) where B is a component of fiamoeba{F). Let ~B* denote the 
negative convex dual of the set B. 

(i) If V ^ -B* then = 0(e^'^l''l) for any (i. 

(a) //x G i? then Or — 0cxp{e^'''^) for all r. 

(Hi) //x G dB but the dot product with r is not maximized over B at x. then Or = Ocxp(e~''"). 

(iv) If r is proper and flcrit{r) is empty, then a^ = Ocxpi^ '^"^)- 

Proof: The first three statements follow directly from the integral formula (|2.ip by taking x • r 
to +00 in (i) and taking x' • r > x + r in (ii) and {Hi). The fourth conclusion is an immediate 
consequence of something we will prove in SectionO under the hypotheses, the contour of integration 
in (|2.ip may be deformed so that Re {— r • y} < — r • x for every y on the contour. □ 

2.6 Quadratic forms and their duals 

Let S denote the standard Lorentzian quadratic xl ~ X2 — ■ ■ ■ — x'^. Any real quadratic form A with 
signature (1, d) may be written as o AI^^ for some invertible linear map M. We now define the 
dual quadratic form A* in two ways. The classical definition is that ^*(r) is the reciprocal of the 
unique critical value of A on the set r'^-'^^ := {x : r • x = 1}. It is easy to compute the dual S* to 
S. The point x is critical for S\^{i) if and only if VS'Hr, that is, if and only if x||(ri, —r2, . . . , —rd). 
Thus the unique critical point of S\^(i) is (ri, — r2, . . . , ~'r'd)/{r1 — r"^ — ■ ■ ■ — r'^) and the reciprocal 
of S there is S*{t) :— 'r\ ~ r2 ~ ■ ■ ■ — r^. In other words, S* in the dual basis {ri, . . . ,rd} looks 
exactly like 5* in the original basis {xi, . . . , Xd}- For the second definition, note any real quadratic 
form A with signature (1, d) may be written as So for some invertible real linear map M . Let 
M* denote the adjoint to M, that is, < A/*r, x >=< r, Af x >; in our coordinates, this is just the 
transpose. We see from the diagram below that Mx is a critical point for A\^{i) if and only if x is a 
critical point for S'|(s*r)'j leading to the alternative definition A*{y) = S*{M*r). 

For computation, it is helpful to compute the matrix for the quadratic form A*. We have 

A(x) = S{M-^x) = y.^ [M-^f DM-^y. 
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where D is the diagonal matrix with entries (1, — 1, . . . , — 1). Thus the matrix for A is (Af ^Y" DM ^ . 
On the other hand, since A*{r) = S'*(M^r) = vMDM^v'^, we see that the matrix for A* is MDM"^. 
In other words, the matrices for the quadratic forms A and A*, are inverse to each other. 

Our definition of the dual quadratic is coordinate free in the following sense. Let A = SoM~-^ as 
above, and let v — (wi, . . . , Vd) denote coordinates in which A is represented by the standard form; 
in other words, v = M^^x and 

A^S{M-\^))^vl-vl vl 

Suppose that an element L e (M'')* is represented by {£i, . . . ,£d) in v-coordinates, that is, L 
maps '^a.jVj to '^aj£j. Then Lx = i, . . . £rf)Af ^^x, that is, L is represented by the row vec- 
tor {ii, . . . ,£d)M~^ with respect to the x-basis. Computing in the x-basis, using this row vector for 
L and the representation MDM^ for A computed above, we have 

A*{L,L) - {£i,...ed)M~^{MDM^){M-Y{h,---,edf 
= {£i,...£d)D{£i,...£df . 

In the V coordinates, A — S and A* — S*, whence A*{L,L) — £\ — £\ — ■ ■ ■ — £'^ and we see that 
dualization indeed commutes with linear coordinate changes. 

Dual quadratics are important because they and their partial derivatives appear in the asymptotic 
formulae for given in Theorems 13. 7[ 13.91 and 16.91 In order to interpret such asymptotic estimates 
and series, it is good to know the size of A* and its partial derivatives. It is easy to see that if F 
is homogeneous of degree n then dF/drj is homogeneous of degree n — 1. It follows that for any 
multi- index m e (2+)"^, the m-partial derivative of (A*)" is homogeneous of degree 2a — |m| and 
hence that 




(2.12) 



The upper estimate is sharp, in the sense that the left-hand side is 0(|rp" I™') except on a subset 
of positive codimension where the m-partial derivative may vanish. 

2.7 Linearizations 

The Fourier integral in (|2.ip turns out to be much easier to evaluate if the function / in the de- 
nominator is replaced by its leading homogenous part. Unfortunately, if (7 is a polynomial with 
homogenous part q, then the fact that g — (7 is of smaller order at the origin than q does not imply 
that q ^ q, which would be necessary for a straightforward estimate of q~^ by q~^ . However, on 
any cone where q does not vanish, we do have such an estimate, and in fact a complete asymptotic 
expansion of q~'' in decreasing powers of q. 
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Lemma 2.24 (expansion in decreasing powers of one function). Suppose that q{x.) = q{x.) + i?(x), 
where q is homogeneous of degree h, and R is analytic in a neighborhood of the origin with i?(x) = 
0(|x|''+^). Let K he any closed cone on which q does not vanish. Then on some neighborhood of 
the origin in K , q does not vanish and there is an expansion 



n=0 



c(m,n)x" 

.\>n{h+l) 



(2.13) 



0(|x| 



(2.14) 



Furthermore, 

g(x)-^- c(m,n)x™q(x) — 

\in\-hn<N 

on K as x—f 0. An expansion of the same type is possible /or p(x)g(x)^'* whenever p is analytic in 
a neighborhood of the origin. 



Proof: Let i?(x) = J2\m\>h+i b{m)x"^ be a power series for R absolutely convergent in some ball 
centered at the origin. Let 

sup|,|g5^EIKm)||x|'" 



inf|x|eas,nK g(x) 
Then by homogeneity, 

^ |b(m)x-"| 
^ |g(x)| 

on the e/(2A/) ball. The binomial expansion (1 + u)^'' = X)n>o Cn')'^^ converges for |u| < 1 and in 
particular for \U\ — 1/2. Therefore, plugging in fe(m)x™/g(x) in for u yields a series 



ri>0 



that converges on i?e/(2M) ri K. Multiply through by g ''to get (|2.13[) . Convergence on any neigh- 
borhood of the origin implies the estimate ()2.14p . □ 



3 Results 

3.1 Cone point hypotheses and prehminary results 

We are interested in the asymptotics of the power series coefficients of a rational generating 
function Fq, in cases where there is a cone singularity and previous known results do not apply. 
Among the properties of Fq discussed in Section [2] there are a number of hypotheses and notations 
that will arise repeatedly. So as to be able to refer to these en masse, we state them here. 
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Hypotheses 3.1 (quadratic point hypotheses). 



1. Let F be the product PqF^^ ■ ■ ■ Fp^ of an analytic function Pg with nonzero real powers of 
Laurent polynomials Fj with no common factor. Assume without loss of generality that Sj ^ Z"*" 
(since otherwise we may absorb Fj^ into Pq). 

2. Let B be a component of the complement of flamoeba([Yj^-^ Fj) so that F has a Laurent series 
expansion on B. 

3. Let Y be a dual vector in the dual cone —B* and assume r is proper with r • x maximized at 



4- Assume that fiW{r) is finite and nonempty. Let w G Tr be an element of flW{r) and denote 

z := Xmax + iw , Z := exp(z) . 

The remaining assumptions enforce a particular set of degrees for the denominator, namely a 
real power of a quadratic together with positive integer powers of smooth divisors. 

5. With Z fixed, we let P be the product of Pq with all Fj such that Fj{Z) ^ and collect terms, 
writing 

F 



Denote q := Q o exp, hj :— Hj o exp, p := P o exp, and denote the homogeneous parts of q and 
hj by q := fihom{q, z) and hj := phom{hj, z). 

6. Assume that q is an irreducible quadratic with signature (1, —1, ... — 1) and let M be a linear 
map such that q — S o AI^^. We allow s — 0, in which case there is no quadratic factor 
vanishing at z. 

7. Assume that hj are linear and that Uj are a positive integers. 

Remark. We lose Uttle generality in assuming fiW(r) is non-empty in clause 4 above, for the following 
reason. If fiW(r) is empty, then part {iv) of Proposition 12.231 guarantees that |ar| is less than x"'' 
by a factor that grows exponentially with |r|. 

In the Aztec and cube grove examples, at the point Z of interest, s = 1, fc = 1, rii = 1, in other 
words, the denominator of P is a product of (the first power of) a quadratic and a smooth factor. 
In the QRW example, rj = 2 but fc = 1 (at each of the two quadratic points, only one of the other 
factors vanishes). There are contributions at the quadratic points (where s — k — ni — 1) but 
they turn out to be dominated by the contributions at smooth points (s = 0). In the superballot 
example, rj = 2 with Fi ~ \ — Axz, F2 = 1 — x — y — z + Axyz, si — —1/2 and S2 = —1. At the 
quadratic point, Z — (1/2, 1/2, 1/2), F2 is quadratic, s = 1, and ni = 1/2. In the graph polynomial 
example, 77 = 1 and s = (3. 
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We extend the expansion in Lemma 12.241 to the generahty of the quadratic point hypotheses as 
foUows. 

Lemma 3.2 (general quadratic point expansion). Assume the quadratic point hypotheses. Let K 
he any closed cone on which g Jlj=i '■^ nonvanishing. Then there is some neighborhood of in K 
such that for all N > 1 the following estimate holds uniformly: 

k 

/(x,„a. + zw + y)- c(m,^,n)y™q(y)-^-^I]/i,(y)-"^-" + 0(|yp^+H+^) . 

m,t,n:\m\-2i-kn<N 3 = 1 

(3.1) 

The sum is finite because c(|m|,^, n) vanishes unless |m| > 3£ + (fc + l)n. 



Proof: Apply Lemma [2. 241 once with g(x + •) in place of q, for the given value of s, and once with 
n^=i ft-j'^(x+ •) in place of q, setting s = 1. This yields two convergent power series. Multiply the 
two series together and multiply as well by the power series for p(x + •). □ 

The results in this paper can be summarized as follows. First, a-r is well approximated by a 
sum of contributions indexed by 6W(r), these contributions being integrals localized near points 
Xmax + iw, for w G fiW(r). Secondly, depending on the geometry at Xmax + *w, this contribution is 
well approximated by a certain explicit funtion of r. The result giving the decomposition as a sum is 
stated as Theorem 13 . 31 below . with the remaining theorems in this section giving the contributions in 
various special cases. It should be noted that Theorem 13.31 is like a trade for the proverbial "player 
to be named later" , in that it allows us to state a complete set of results even though the meaning 
will not be clear until the other results have been stated. 

Theorem 3.3 (localization). Assume the quadratic point hypotheses and notations. Then there is 
a conical neighborhood M of r in {W^)* and there are chains {Cw : w G fiWijQ)} defined in the text 
surrounding Theorem \5.4\ in Section below, such that 

Or = ^ ficontrib{w) + Ooxp (x^L) ■ (3-2) 
v,ef!W{ro) 

The estimate is uniform when |r| cx3 while remaining within M. The summand is defined by 



f^contr^bM := { — / ^.^ ' dz . (3.3) 

j^^/ Jc„ '?(z)'' n,=i ^^A^)^ 



Proof: This is an immediate consequence of Corollarv 15.51 below. □ 



3.2 Asymptotic contributions from quadratic points 

Next, we identify the contributions 6contrib(w). In the case where w is a smooth point (s = 
0, fc = l,ni = 1), these are already known. A formula involving the Hessian determinant for a 
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paramctrization of the singular v ariety V o f F was given in PW02I . Theorem 3.5], which was then 
given in more canonical terms in BBBPO 



PW02 



BBBPOSj). Assume the quadratic point hypotheses and suppose that 



Theorem 3.4 ( 

0, fc = 1 and ni — 1, so 7i is a simple pole for F. Let 

denote the gradient of H := Hi in logarithmic coordinates and let k — k{z) denote the (possibly 
complex) Gaussian curvature ofVf at z. Suppose that n ^ 0. Letting \ ■ \ denote the euclidean norm, 
we have: 



flcontnb{w) - (27r|r|)^^"''^/^ 



p{z) 



vMz)|Viog| 

The estimate holds uniformly over a sufficiently small neighborhood of r such that: (i) the quadratic 
point hypotheses are satisfied, (ii) k ^ 0, and (iii) the point Z = Z(r) varies smoothly. The square 
root should be taken as the product of the principal square roots of the eigenvalues of the Gauss map. 
□ 



In the case where w is on the transverse intersection of smooth (local) divisors, formulae are 
also already known. There are a number of special cases, depending on the dimension of the space, 
the dimension of the intersection, and the number of intersecting divisors. We will not need these 
res ults in this paper (we need only the upper boun d in L emma 15. 9|) but statements may be found 



in [PW04| . Theorems 3.1, 3.3, 3.6, 3.9, 3.11] and in jBP04| . The novel results in this paper concern 
the case at a quadratic point, that is, where s ^ 0. Let N*x(i?) denote the dual to the tangent cone 
tanx(-B). The cone N*x(S) will have nonempty interior. By contrast, in example 12.191 the cone 
tanx(-B) is always a half space and N*x(i?) is always a single ray. 

Definition 3.5 (obstruction). Assume the quadratic point hypotheses and notations. Say that r is 
non- obstructed if r is in the interior o/N*x(i?) and if for any x in the boundary of the cone of 
hyperbolicity of A, the cone K'^'^(x) contains a vector v with r • v > 0. 



This condition is not transparent, so we pause to discuss it. First, note that the non-obstruction 
condition will turn out to be satisfied for all r in the interior of N*x(S) when fc = (locally, the 
denominator of F is an irreducible quadratic). To see this, recall from Proposition 12 . 71 that the cone 
of hyperbolicity of g is a component of its cone of positivity. At any point v on the boundary of this 
cone, other than the origin, q is smooth and hence fihom(g, v) is linear, vanishing on the tangent 
plane at v to {q = 0}. The normals to these planes are precisely the extreme points of the cone 
N*x(i?). Therefore, for any r in the interior of N*x(i3), r is not perpendicular to the tangent plane 
at to {q = 0} at any point v ^ 0, which implies that r is non-obstructed. An example where there 
are obstructed directions interior to N*x is as follows. 
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Example 3.6 (obstruction). Suppose the denominator of F is H := H1H2H3 := (1 — — y)(l — 
XY). Then h := Cxy{x + y). The cone tan(Q g)(i?) is the negative orthant. The dual cone N*x(-B) 
is the positive orthant. The vector r — (1, 1) lies in the interior of the dual cone. Let x = {t, —t) for 
some t ^ 0. Then K'^'^(x) is the halfspace {{x,y) : x + y < 0} and r • {x,y) is maximized at zero 
on this cone. 

Secondly, we see that the condition of non-obstruction is not merely technical, but is necessary 
for the conclusions we wish to draw. To elaborate, we would like our asyniptotics to be uniform 
as r varies over the interior of N*x(-B)- Unfortunately, this is not always possible. In the previous 
example, if F — 1/H = 1/[(1 — ^)(1 — — XY)], then ~ min{ri,r2}. Analytic expressions 
for Or will not be uniform as r approaches the diagonal. This is in fact because movement of the 
contour of integration in (|2.ip will be obstructed, requiring different deformations for r in the positive 
quadrant on different sides of the diagonal. 

Theorem 3.7 (quadratic, no other factors). Assume the quadratic point hvvotheses \3.1[ and suppose 
that k — 0, in other words, F = P/Q'^ with no further factors in the denominator. Assume further 
that 

s 7^ 0, -1, -2, . . . and s^ --\,--2,... . 

Let c(m, n) be the coefficient o/x™(j(x)^^^" in the expansion (|2.13p . Let K* he any compact subcone 
of the interior oflSI*. Then, uniformly over r e K*, there is an expansion 

I /If I am . , 

ficontnb{w) ^ ^ ^ ^ ' ttZ^'V V c(m, n)(-l)l'"l ^ ( r(r)"+"^''/' ) . 

y ^ y I ' n |ml>3n 

(3.4) 

The series is asymptotic in the following sense. For any N , restricting the series to terms with 
|m| ~ 2n < N yields an approximation whose remainder term is 0{\r\'^'^~'^~^), all of whose terms 
are generically of order |i-|2s-<i-A^+i_ /Jp(Z) then 

Remark. Comparing to equation 12.121 we see that the remainder terms are no larger than the first 
omitted term of (I2.12p . For a true asymptotic expansion, this should be smaller than the last term 
that was not omitted, but in general there may be directions r in which (9™/9r™)g*(r)*~''/2 jg of 
smaller order than |r|2s-'*-|i"l . This may occur after the first term in the expansion p.4p . though 
not in the leading term (|3.5p . Also, by Theorem l3.3[ we may be adding up several of these formulae, 
thereby obtaining some cancellation. For example in the case of the Aztec diamond, — when 
^ r j is odd. This manifests itself in the symmetry F{Z) = F{—Z), and in two quadratic points at 
(1, 1, 1) and (—1, —1, — !)• Contributions from the two quadratic points will sum or cancel according 
to the parity of r. 

As a corollary, for ease of application, we state the asymptotics in the three variable case for a 
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single power of Q in the denominator. Theorem 13.71 is proved in Section 16.41 while Corollary 13.81 
follows immediately. 



Corollary 3.8. Assume the quadratic point hypotheses with d ^ 3, k ^ and s — 1. Let c(m, n) 
he the coefficients in the expansion \2.1SX Let K* be any compact subcone of the interior of N* , the 
dual cone to tanx(i?). Then, uniformly over r S K*, there is an expansion 



ficontnbiw) - LUz-'-^ ^ c(m,n) 

n=0 |m|=Ti 



q*(r) 



-1/2 



(3.6) 



Here, asymptotic development means that if one stops at the term n ^ N — 1, the remainder term 
will be Odrl^""^^^), while the last term of the summation will be of order |r|^^. Ln particular, if 
P(Z) ^ then 



f^contr^bi^) ^ ^iMz- 
27r 



r(r) 



-1/2 



uniformly on K* 



(3.7) 

□ 



Remark. Again, the leading term estimate p.7p is a true asymptotic estimate, while the right-hand 
side of p.6p may vanish for certain m and r. 



3.3 The special case of a cone and a plane 

Our last main result addresses the simplest case where the are both a quadratic and a linear factor. 
The case of a quadratic along with multiple linear factors is also interesting. We address this in 
Section [^31 Because there are a great number of subcases and we have no motivating examples, 
we do not state here any theorems about that case, and instead describe in Section 16.51 a number 
of results that pertain to this case. In the case of a single factor of each type, in three variables, 
significant simplification of the general computation is possible. The remaining results concern this 
special case. 

Assume the cone-point hypotheses with = 3,s = 1 and k = 1. Because fc = 1, we drop the 
subscript and denote H :— Hi. We assume also that the linear factor i := hi of the homogeneous 
part of (QH) o exp shares two real, distinct projecive zeros with the quadratic factor and we 
denote these by ai and a2- The given component B on which the Laurent series arZ"" converges 
is the intersection Bi O B2 oi a. component of 6amoeba(Q)^ and 6amoeba(i?)^. By hyperbolicity, we 
know that the quadratic g is a scalar multiple of a real hyperbolic quadratic; multiplying by —1 if 
necessary, we may assume the signature to be (1, 2); in particular, we may write 

q(Av + w) = X"^ - |wp 

for some v E tanx(i?i) and all w G tanx(i?i)^. The set Bi is a cone over an ellipse £ and its dual 
~Bl is a component of positivity of the cone q* = 0. The linear function h may be viewed as a 
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point of (R^)*. Figure [5] shows a plot of q* = and of the point h in (R^)*. Also shown is the hne 
of points r for which q*{r, h) = 0. These shapes in the projective (r|s|i)-space (RP^)* are shown via 
their slices at i = 1. 




Figure 8: the cone N* depicted by its slice at t = 1 



The assumption that ai,a2 are real implies that the point h lies outside B^. The normal cone 
N*x(S) is the convex hull of the normal cone of Q and the normal cone {h} of H. This teardrop- 
shaped is the entire shape shown in Figure [8l The tangent lines to B* from h intersect B* in two 
projective points, namely those r for which q*{r,h) = q*{r,r) = 0. The non-obstructed set is a 
disjoint union i3* U E, where the cone E is the non-convex region N*x(i3) \ Observe that the 
dotted arc in figure [8] is obstructed and thus is in neither B^ nor E, these being the two components 
of the non-obstructed set. 

To state the final theorem, we must define one more quantity. If A is a homogeneous quadratic 
and L is a linear function, define a quantity Res'^' as follows. Let denote the form (z dxdy — 
y dzdx + X dydz) /{A-L). The second iterated residue of is a 0-form, defined on the two lines ai , 02 
where A = L = 0. Because 9 is homogeneous of degree zero, its second residue has a common value 

at any affine point in the line Uj. We let Res^^^ = Res^''^(aj) denote this value. In coordinates, we 

(2) 

have a number of formulae for Res^ , one being 



Res(^)(a) 



dA dL dA dL 

dx dy dy dx 



(3.8) 



{x,y,z)<^c 



Theorem 3.9 (quadratic and one smooth factor). Assume the quadratic point hypotheses with 
d = 3, s = 1 and k = I and let £ denote the linear factor, hi at the point z. Assume p[z) ^ 
and assume that the two projective solutions ai,a2 to £ = q ^ are real and distinct, so that the 
non-obstructed set N* is the union of an elliptic cone B^ and a nonconvex cone E as describd above. 

Let q* denote the dual to the quadratic q. Let arctan denote the branch of the arctangent mapping 
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(0,oo) to (0,7r/2), while mapping {~oo,0) to (7r/2,7r) rather than to (— 7r/2,0). Then 



j]contrib{w) = Z-^PiZ) 



■ arctan 



R 



(3.9) 



where the remainder term satisfies R = 0{\r\ ^) uniformly as r ranges over compact subcones of 
Bl. On the other hand, we have the estimate 

ficontrib^ = Res(^)p(Z)Z"'' + R 

where R = 0(|r|~"'^) uniformly as r ranges over compact subcones of E. 



4 Five motivating applications 



One feature is common to all but one of our applications, namely that is on the boundary of the 
amoeba of the denominator of the generating function. In this case, by part (Hi) of Proposition l2 . 23l 
the coefficients Or decay exponentially as |r| — s- oo in directions f for which supyg^j^j^^^^-) f • y > 0, 
in other words for r ^ N*, the dual cone to tano(i?). In such a case, the only significant (not 
exponentially decaying) asymptotics are in directions in the dual cone (tano(i?))*. We therefore 
restrict our attention in every case but the superballot example to r e (tano(S))*, and consequently, 

to X^iax 0. 



4.1 Tilings of the Aztec diamond 



The model 



The Aztec diamond of order n is a union of lattice squares in Z^, whose vertices are the pairs (±i, ±j) 
with i, J > 1 and i + j — n or n + 1. Thus the order 1 Aztec diamond consists of the four squares 
adjacent to the origin, the order 2 diamond consists of these together with the square centered at 
(3/2, 1/2) and its se ven images under the symmetries of the lattice rooted at the origin. This was 
defined in EKLP92l | , where questions were considered regarding the statistical ensemble of domino 



tilings of the Aztec diamonds. A domino tiling of a union of lattice squares is a representation of 
the region as the union of 1 x 2 or 2 x 1 lattice rectangles with disjoint interiors. Figure [9] shows an 
example of a domino tiling of an or der 4 Azte c diamond. The set of domino tilings of the order n 
Aztec diamond has cardinality 2© |EKLP92t . Let fin be the uniform measure on this set, that is. 



the probability measure giving each tiling a probability of 2 -(3) . Limit theorems for characteristics 
of /i„ have been proved, the most notable of which is the Arctic Circle Theorem which states that 
outside a (1 + e) enlargement of the inscribed circle the orientations of the dominos are converging 
in probability to a deterministic brick wall pattern, while inside a (1 — e) reduction of the inscribed 
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Figure 9: the Aztec diamond of order 4, tiled by dominoes 



circle the measure has positive entropy 
locations inside the circle were given in 



JPS9S L A new proof and a distributional limit at rescaled 



CEP9a Theorem 1]. 



Via an algorithm called domino shuffling ProOSj, the following generating function was found. 
Color the square centered at (« — ^, j — ^) in the Aztec diamond of order n black if i + j ' + n is 
odd and white if i + j + n is even. A domino is said to be northgoing if its white square is the 
(0, l)-translate of its black square. For i + j + n odd, let p{i,j, n) denote the probability that the 
domino covering the square centered at (i — ^, j — ^) is northgoing. Then the generating function 
for the quantities p{i,j,n) is given by (jl.3p . which we recall here: 



Z/2 



(1 - YZ)[1 -{X + x-i + Y + y-i)z/2 + ^2] 



The sum is taken over n > 1 and —n < i, j < n with |i — ^| + |j — ^| < n and i + j + n — 1 = 
modulo 2. The first results on these probabilities were derived using bijections and ot her algebraic 
combinatorial methods jEKLP92i] . The formula (fT31) for F was first proved in [jPS98| . 

Related to the probabilities p{i,j,n) are the creation rates c{i,j,n) := 2{p{i,j,n) — p{i,j — 
1, n — 1)). Their generating function is the same except that the factor (1 — YZ) in the denominator 
is omitted: 



create 



Z/2 



1-{X + x-i + Y + y-i)z/2 + Z2 



CEPga, this 



The creation rates have an interpretation via the shuffling algorithm [GIP96|. In 
generating function is related to Krawtchouk polynomials, and the Arctic Circle limit theorem 
is derived by combining some known fact s about these polynomials with contour integrals that 
determine coefficients of these polynomials CEP96| . The method of contour integration used there 
is in some sense a special case of the contour methods in the present paper, tailored to the particular 
case in which one need compute coefficients only of the relatively simple expression {1 — Z)'^~^ {1+ Z)'' . 

It was guessed that the asymptotics of the crea .tion rates c(i, j, n) could be derived directly 



from the generating function F along the methods of PW02I : |PW04| , and some progress along these 
lines, was made by Cohn and Pemantle in an unpublished manuscript. Via resolution of singularities. 
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they reduce the Cauchy integral for J^create which the denominator factors into irreducible, 

smooth divisors. The phase function in (|1.8p becomes e~'^^^^ for some (f> that is more complicated 
that r • X. Technical difficulties with the phase function and the topology of the integral prevented 
Cohn and Pemantle from completing their proof. 



The amoeba and normal cone 

We apply the results of Section [3l An outline is as follows. After verifying the quadratic point 
hypotheses, the localization Theorem 13.31 computes asymptotically as a finite sum 

ficontrib(w) . 

wefiW(r) 

The point (0, 0, 0) is on the boundary of the component B and is in fact a quadratic point. We will 
compute its normal cone N* which is the teardrop shaped region shown in figure [231 Outside of N*, 
the probabilities decay exponentially. When r £ dN* we cannot say anything, but for r interior to 
N* we will obtain, via Theorem 13. 9[ a 2-periodic contribution at the critical points ±(1, 1, 1). The 
leading term asymptotics (|4.2p will follow once we show all other contributions to be negligible. 



Corresponding to the notation in the quadratic point hypotheses, we write F = P/ (QH) where 
Q :=1-{X + + Y + Y-^)Z/2 + Z^, H ■=!- YZ and P := Z/2. Using a computer algebra 
system to compute a Grobner basis for {Q,QxtQytQz}: we find that Vq is singular precisely 
at Z = ±(1,1,1). Letting q := Q o exp and q := hom(q, 0), wc find at the point (1,1,1) that 
q{x, y, z) — — ^x^ — iy^; the computations for the point (—1, —1, —1) are analogous and are done 
at the end of the discussion. We see that near (1, 1, 1), Q is an irreducible quadratic, while h is linear, 
with linearization h{x, y,z) = y + z. To specify B, observe that the components of fiamoeba(/)'^ are 
intersections of complements of fiamoeba((5) with components of the complement of 6amoeba(7J). 
A glance at the series J^Pihjin)^^^''^"' shows that the series is convergent for any fixed X and 
Y as long as Z is sufficiently small. Hence the component B of the complement of fiamoeba(Q7f) 
corresponding to this series is the one containing (0,0, —A) for sufficiently large A. The amoeba of 
1 — YZ is just the line y — —z in log space, and the component of j3amoe ba(iif ) co ntaining the ray 



(0, 0, —A) is the halfspace Bi {y + z < 0}. Turning to Q, we recall from [GKZ94 Chapter 6] that 
the components of the complement of the 6amoeba((5) correspond to vertices of the Newton polytope 
SP{Q). The ncwton polytope is an octahedron with vertices (±1, 0, 1), (0, ±1, 1) and and (0,0, 1±1). 
There is one vertex, namely (0,0,0), for which (0,0,— A) is in the interior of the normal cone. Let 
B2 to be the component of 6amoeba((5)'^ containing a translate of this cone. Let B — BiC] B2- This 
completes (l)-(2) of the quadratic point hypotheses. 

As discussed at the beginning of Section 21 in the case where S dB, we will be chiefly interested 
in asymptotics in directions r for which r -x < for x G B. Let us verify that G dB. Let Z — U XY. 
Observe that if X,Y,U < I then the series p.3p is absolutely convergent. Sending U, X,Y tol sends 
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{X, Y, Z) to (1,1,1) which is therefore on the boundary of the domain of convergence of (|1.3p : hence 
= log(l, 1, 1) is on the boundary of fiamoeba(Qi? ) . We now compute N* := —{K'^'^'^)*. This was 
done in general in Section 13.31 so to complete the description, we need merely to identify the dual 
quadratic q* and the dual projective pont h. The quadratic q is already diagonal: q — z^ — (a;^+y^)/2; 
hence q* — {l/2)t^ — — . Letting (f , fi,") be the unit vector r/|r|, we obtain the plot in figure [TUl 
The projective point h is the point (0|1|1), which in the t = \ slice is given by (0, 1); this is outside 



h = (0,1) 




Figure 10: the disk B2 and the point h 
the dual cone B^, reflecting the fact that q and h have two common real solutions. 



Classifying critical points 

Suppose that r is in the interior of i?2- Then x,nax — and r is non-obstructed and proper. To finish 
verifying quadratic point hypotheses (4)-(5), we need to identify 6W(r) and check that it is finite. 
As noted before, it will turn out that X)wG6W(r) ficontrib(w) is dominated by the contributions 
from w = and w ~ {tt,tt,tt). We may therefore identify the remaining critical points somewhat 
less explicitly. 

Finding the critical points requires an explicit stratification of Vp- The coarsest Whitney strat- 
ification is as follows. 

Vi := {(1,1,1), (-1,-1,-1)} 
V2 := VqHVhXVi 
V3 :- Vh\(ViUV2) 
V4 Vq\(ViUV2) 

defines a Whitney stratification of Vf- The points of Vi are isolated (quadratic) singularities of 
Vq, while the remaining strata are Vq, Vh and their interstecion, which may be parametrized by 
{(z^^, z^^ , z) : z £ C*}. By definition, any function is critical on a zero-dimensional stratum, whence 
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both points of Vi are critical for all r G N*. Below, we will show that in fact ficontrib(w) = 6(1) 
for exp(iw) e Vi. When r is in the interior of N*, we will show that the remaining critical points 
break down as follows. 

V2 : No critical points 

V3 : No critical points (4-1) 

V4 : Finitely many critical points 
By Theorem l3.4| the critical points in V4, which are smooth, each contribute o(l) to the asymptotics, 
so we will be done once we evaluate the contributions from ±(1,1,1) and prove (|4.ip . 

Turning to the issue of counting critical points, we begin with the easiest stratum V3. Recall 
from (|2.10p that on the smooth stratum Vh, the point Z is critical if and only if ViogZ is parallel 
to r. The logarithmic gradient of H the constant vector (0, 1, 1), which is on the boundary of N*, 
whence V3 contains no critical points interior to N*. To compute critical points on V2, we evaluate 
Viog(3(z*^, z~^, z) and find that independent of z, we always obtain the projective point (t^j \^ !)■ 
This shows that Vq intersects Vh transversely, and by (|2.1ip . that V2 produces critical points only 
when r is in the union of two projective lines, one joining (0, 1) to (^, ^) and the other joining (0, 1) 
to (—1/2, 1/2). This union does not intersect the interior of N*. 

To solve for critical points in V4, fix r = (r, s, t) and solve the equations (|2.10p : Q = 0, tXQx — 
rZQz = and tYQy — sZQz = 0. Multiplying each of these by 2xy clears denominators and 
allows us to use a computer algebra system to compute a Grobner basis for the solution. With 
lexicographic term order plex(a;, y, 2), the almost-elimination polynomial for z is yz(l + z)^(l — z)^ 
times a quadratic polynomial in r, s, t and z^: 

{r* -2r^s^-2 t^r^ + - 2fs^ + t^) + {-2 + 4r\s^ - 2r^ - A th^ - 4 + 2 t") 
+ (r* - 2 r^s^ - 2 t^r^ + 8^-2 fs^ + t*) z* 

It is easy to check that for every r, s, t, this polynomial is not identically zero, hence there are only 
finitely many solutions. The basis contains a polynomial in y and z (over C{r,s,t)) that is linear 
and non-constant in y, implying that for each z there is at most one y. The same is true for x if we 
use the term order plex(j/, x, z). It follows that there are finitely many critical points in V4 for each 
r. Summing up, have verified (4)-(5) of the quadratic point hypotheses for r interior to N*. 

Computing the estimate 

The computation for w = (tt, tt, tt) (hence Z = (—1, —1, —1)) are almost identical to those for w = 
and Z = (1, 1, 1). We do the latter computation and indicate changes needed to do the former. We 
observe also that Q and H are invariant under {X,Y,Z) i—f {—X,—Y,—Z), while the numerator, 
Z/2, is odd; this corresponds to the parity constraint of p{i, j, n) vanishing when i + j + n is even. 

The quadratic point hypotheses have now been spelled out and verified. Let w := and Z := 
(1, 1, 1). To check that we are in the case covered by Theorem 13.91 we need to check that the two 
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projective solutions to q = h = are real and distinct. This is easy: plugging in y = — z, we get 

1 2 _ 1, 
2^ 2' 



— — = which has the two real solutions y = —z ~ ±x. 



The quantity q*{r,r) is the quadratic that is positive on the interior of the disk, reaching a 
maximum of 1 at (0, 0, 1) and vanishing on the boundary of the disk. In coordinates, it is given by 

g*(r,r) = f - 2s'^ - 2r^ . 

The quantity q*{r, h) is equal to t — 2s. This vanishes on the line shown in figure [3731 The branch of 
the arctangent chosen in the conclusion of Theorem l3 . 91 varies continuously through 7r/2 as r, h) 
varies through zero and the argument of the arctangent passes through ±oo. The arctangent goes 
to zero where g*(r, r) = and g*(r, ft,) > (the part of the boundary of the disk to the left of the 
vertical line) and to tt where g*(r,r) = and q*{r,h) < (the part of the boundary of the disk 
to the right of the line). The residue Res'-^-' is immediately computed from (|3.8p and is equal to 1. 
Finally, we have P{z) = 1/2 and q*{h, h) = —1. Thus, as r varies over the interior of the projective 
disk B2 we have 



1 Vq*{r, r)\ 1 / Vt^ - 2r^ - 2s2 \ 

DContrib(O) ~ — arctan — = — — — arctan 

^ ' 2t: y q*{r,h) J 2Tr ^ t ~ 2s j 

The computation for w = (7r,7r,7r) is entirely analogous, leading to the same contribution but with 
an extra sign factor of (— l)*+J+"+i. We have already shown that all other contributions are of order 
0(|r|~"'^). Therefore, we may sum these results to obtain the formula 

(r + s + i)mod2 f - 2r'^ - 2s'^\ , , 

- - ' arctan (4.2) 



1 — 2f2 — 2s^ \ 

which is equal to tt^^ arctan I on odd multi- indices, and zero on even multi-indices, 



whe re (f, s] := (r/t,s/t). This formula, appearing as [CEP961 Theorem 1], is plotted in figure [TTl 



see 



CEP961 Figure 2] for a contour plot of the same function. 



4.2 Cube groves 
The model 



CS04I : IPS04 

we define a collection of lattice subgraphs known as cube groves. Let L„ be 



After 

the triangular lattice of order n > 0, by which we mean the set of all triples of nonnegative integers 
(r, s, t) e (2+)'^ such that r + s + t — n with edges between nearest neighbors (thus the degree of an 
interior vertex is 6). We depict this in the plane as a triangle with n + 1 vertices in the top (zeroth) 
row, and so on down to 1 vertex in the 71*'* row. 

The cube groves of order n are a subset C„ of the subgraphs of L„. The set C„ has a description 
where one begins with the unique cube grove of order zero, then produces sequentially groves of 
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-0,7 



Figure 11: asymptotics for northgoing probabilities in the Aztec Diamond 
orders 1, 2, . . . , n, each produced from the previous by a "sh uffle" which injects some information 



in a manner similar to the domino shuffling used by CEP96l | in studying and enumerating domino 
tilings of the Aztec diamond. The set C„ has other, static definitions in terms of graphs that look 
like stacks of cubes and in terms of graphical realization of certain terms of generating functions 



(see PSOSj ) , but here we will take the shuffling procedure to define the set C„ of order n cube groves. 



Define Cq to be the singleton whose element is the one-point graph. If T is a downward-pointing 
triangular face of L„, let T' be the rotation of T by 180° about its center. The union of the vertices 
of the triangles T' is a translation of the graph Ln+i, provided that one adds in the three corner 
vertices of Ln+i- The edge sets of the triangles T' are disjoint and their union is the edge set of 
Ln+i, provided that one adds in the six edges adjacent to corner vertices. 

• — » — • 



Figure 12: an order 4 cube grove, shuffled to become an order 5 grove 

Given a cube grove G G L„ and a downward-pointing triangular face T of L„, let G{T) be the 
collection of graphs on T' that have: no edges if G has two edges in T; one edge if G has one edge 
e e T, in which case the edge of T' must be the edge of T' parallel to e; two edges if G has no edges 
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in r, in which case any two of the three edges of T' wiU do. Let C{G) be the direct sum of G{T) 
as T varies over downward-pointing triangular faces of L„. That is, choose an element of G{T) for 
each T and take the union of these. Figure [T^ shows an order 4 grove, G and one of the 27 elements 
of C{G). Finally, let Cn+i be the (disjoint) union of the collections C{G) as G runs over L„. 

Looking at a picture of a uniformly chosen random cube grove of order 100, one sees regions of 
order and disorder similar to those of the Aztec Diamond. Let Pn{hj) be the probability that the 




Figure 13: a random cube grove of size 100 



horizontal edge with barycentric coordinates n — i — j) is present in a uniformly chosen cube 
grove of order n. The creation rates E„{i,j) may be defined in terms of the shu ffling p rocedure but 
in this case they satisfy the simple relation En-i{i,j) — \{Pn{i,j) ^ Pn- 
We recall here the explicit generating function (|1.4|) . which is derived in 



[PSOSl . Theorem 2]. 



PSOa . Section 2.2]: 



F{X,Y,Z)^ 



2Z^ 



(1 - Z){3 + 3XYZ - {X + Y + Z + XY + XZ + YZ)) 



2Z2 
HQ 



It is quick to verify that longer factor, Q, in the denominator has a quadratic cone singularity and 
that F therefore is singular on the union of a quadratic cone with a smooth surface. The real part 
of this is pictured in figure [Ml 



The amoeba and the normal cone 

All multi-indices in the generating function are nonnegative, so it is an ordinary generating function 
and B will be the component of 6amoeba(i^)'^ containing the negative orthant. Again, we are chiefly 
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X 



Figure 14: the pole variety of the cube grove generating function 

interested in directions r for which Xmax(r) = 0, these being the directions of non-exponential decay. 
The polynomial Q has a single quadratic point at (1,1,1). To compute tanx„,^^(-B), we intersect 
Bi {(x, J/, z) : z < 0} with the cone B2 of hyperbolicity of Q at (1, 1, 1). Changing to exponential 
coordinates via q :~ Q o exp and computing the leading homogeneous term gives 

(l{x,y,z) = q{x,y,z) + 0{\i,\^) 
where q(x, y, z) := 2xy + 2xz + 2yz . 

It follows that B2 is the cone containing the negative orthant and bounded by {g = 0}. The dual 
quadratic is represented by the matrix 



" 


1 


1 " 


-1 


" -1 


1 


1 








1 








1 





1 


~ 2 


1 


-1 


1 


1 


1 







1 


1 


-1 



hence 

q*{r,s,t) ^ rs + rt + st ~ -(r^ +s^ + i^). 

The dual cone —Bf is the subcone of the positive orthant bounded by q* = 0. Again, the point h, 
which is equal to (0, 0, 1) in the (r, s, t) coordinates, lies outside this cone, and again the solutions to 
q — h — are the solutions to z = = xy which are two distinct projective points, namely the a;-axis 
and the y-axis. We could again depict this by the slice through t = I, viewing B2 as the interior 
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of a parabola in the first quadrant, opening in the Northeast direction and tangent to the axes at 
(1,0) and (0,1), with h at (0,0). It is easier to see what is going on if we change coordinates to 
u := (1, 1, 1)/ \/3, letting denote the complementary space. For r = (r • u)u + r-^ with r-*- € U-^, 
we then have |rp — (r • u)^ + |r^p, whence 

= ,2 ^ ^,2 ^ ,2 _ 1 ^ ^ ^ ,)2 ^ 1 (^2 ^ ^2 _^ ,2) _ 2 ^ 

Thus q*{r,s,t) — when |rp = 3|r^p, or equivalently, |r • up = 2|r^p. Viewing projective space 
via the slice r • u = 1, we see that q* vanishes on the circle centered at the origin of radius \/l/2. 
The projective point h = (0, 0, 1) intersects the slice |r • u| = 1 at (0, 0, -s/S), whose projection to If-^ 
has squared norm 2. This is pictured in figure [TSl In these coordinates, the only difference between 
this figure and that for the Aztec Diamond is that the distance from the origin to the point h is 
twice the radius of the circle, rather than \/2 times the radius, and the tangents subtend an arc of 
120° rather than 90°. 

z-0, x-y 




y-0, x-z 



Figure 15: the dual cone in symmetrized coordinates 



Classifying the critical points 

The stratification is similar to that for the Aztec Diamond generating function. There is just one 
singular point of Q, namely (1, 1, 1). This is on Vh as well. The surfaces Vh and Vq intersect in 
the set {x — z = 1} U {y = z = 1}, which is smooth away from (1, 1, 1), leading to the following 
stratification. 

Vi := {(1,1,1)} 
V2 := VqHVhXVi 
Va Vh\{ViUV2) 
Vi := Vq\(ViUV2) 

The logarithmic gradient of H is parallel to (0,0, 1), which is not in N*, so for r e N* there are 
never any critical points on V3. There are critical points on V4, but we verify as before that there 
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are only finitel many. They are smooth, so by Theorem l3.4[ their contributions are o(l). On V2, the 
logarithmic gradient of Q is parallel to either (1, 0, 1) or (0, 1,1). The logarithmic gradient of H is in 
the t direction, so the span of the two logarithmic gradients is either the r-t plane or the s-t plane. 
Neither of these planes intersects the interior of N* (the planes are tangent to N* at the projective 
points (1,0, 1) and (0, 1, 1) respectively). The for r interior to N*, there are no contributions from 
V2; it remains to compute the contribution from Vi. 



Computing the estimate 

Completing the computation as in the Aztec case, we evaluate Res'-^-' using (|3.8p but switching the 
roles of X and z because only the z-derivative of h is non-vanishing. This gives Res(2) = i. With 
P = 2Z^, and only one contributing point w — (0,0,0), we have q*{r,h) — {r + s — t)/2 and 
q*{h, h) = —1/2, whence Theorem 13.91 gives 



TT 1 (r + s-i)/2 



1 

— arctan 



1 / J2(rs + rt + st) - (r^ + s'^ + m \ 

= — arctan — 

77 \ r + s — t I 

Again, the arctangent is taken in (0, tt) so that as we cross the line t = r + s the arctangent varies 
continuously across 7t/2. 



4.3 Two-dimensional quantum random walk 
The model 

We begin with a brief review on one-dimensional quantum random walk (QRW). In the classical 
simple random walk, the law at time n is a probability measure on Z and the evolution operator 
on this law is (1/2)(t_|_ + (1/2)(t_, where cr+ is the right-shift operator a^fi{n) ~ fi{n — 1) and a_ 
is the left-shift operator a^fi{n) = ^{n + 1). In the quantum world, the law at time n is given by 
the values of |^/'(n)p where the wave function if) is not a positive unit vector in L^(M) but rather a 
unit vector in L^(C). Evolution operators must be unitary. While the shifts a± are unitary, linear 
combinations of these such as (1/2)ct+ -I- (l/2)tT_ are not. 



An idea for constructing a quantum simple random walk, apparently due to Mev96l |. is to enlarge 
the space to := Z x {f/, D}, adding a hidden "spin" variable. To take a step of the random walk, 
first the spin is randomized, then all particles with spin up move one step right and all particles with 
spin down move one step left. A number of choices are available for the operator that executes the 
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randomization of spins. One common choice is the Hadamard coin-flip, B 



The 



V2\l -1 ^ 

terminology reflects the fact that the matrix is a multiple of a matrix with ±1 entries, these being 
known as Hadamard matrices. Under this operator, either state (1,0) or (0,1) becomes an equal 
mix of U and D states. Let A be the operator which maps state (n, U) to {n + 1, U) and (n, D) to 
(n — 1, D). If we begin in state (0, C/), then do S :— Ao B, the particle is in an equal mix of states 
(1, U) and (—1, D). If we measure the position, we will have executed a step of QRW. 

The n-step simple random walk is defined to be the operator 5" :— {ABy\ If the position if 
this is measured at time n, the probability of being at position k is |iS"(fc, C/)p -I- |iS"(fc, D)\'^- Since 
no measurement is made until time n, the various possible coin-flips and moveme nts interfere, both 



positively and negatively, and the result is somewhat complicated. The analyses in [NVOOt lABN+Ol ] 
show that unlike classical simple random walk, QRW spread s out lin early, with location distrubuted 
over the interval [— n/\/2, n/\/2]; see also the review article Kem05 |. 



To define a two-dimensional QRW, we need a four-fold auxilliary state. Denote these four states 
by {N, S, E, W}. Any 4x4 unitary matrix may be used for the quantum coin-flip. The Hadamard 
matrix 

/I -1 -1 -1 \ 

-1 1 -1 -1 
-1 -1 1 -1 

V -1 -1 -1 1 / 



-4 



is known ( http: / / www.santafe.edu^moore / gallery. ht]XLl | ) as the Hadamard quantum coin-flip. N.B.: 
This is different from the Hadamard QRW in [GJS04| , which also uses a Hadamard matrix, namely 
the tensor product of two copies of the one-dimensional Hadamard matrix. A step of the two- 
dimensional Hadamard QRW is the product AU where A maps ((r, s), N) to ((r, s + 1), N), and so 
forth. The following generating function for the probab ility amplitudes of a QRW in any dimension 
with any quantum coin-flip matrix is given in BPOTl Proposition 3.1]. Let M be obtained from 
U by multiplying the first row by XZ, the second by YZ, the third by X~^Z and the fourth by 
Y^^Z. We consider the rows and columns of M as indexed by the ordered quadruple {E, N, W, S). 
Then an entry of M* such as {M'^)m,e counts the number of t-step paths from N to i?, weighted 
by M: the Z* coefficient of this is the wave function at position ((r, s), N) and time t starting 

from ((0,0), iV). Summing in t shows that the components of (/ — M)^^ = ^ Af (r, s, t)X''F*Z* 
are the generating functions for the wave function at all positions and times: each M(r, s, t) is a 
matrix, whose (^,r7)-entry is the generating function ^ ^ c(^, 7y; r, s, i)X'"y*Z* where c(r, s,i) is 
the probability amplitude, starting from state ((0,0),^) at time zero, of being in state ((r,s), 77) at 
time t. 

The entries of (/ — M)^^ have denominator (1 — Z^)Q where 



:= 1 - 2- 



x + x- 



Y + Y- 



-Z + Z^ 
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We recognize the same polynomial factor that occurred in the Aztec denominator. The numerators 
in the first row are half of the following. 

Pi = 2- {Y + Y-^ + X-^)Z + 

P2 = XZ - {1 + Y-^X)Z^ + Y-^Z^ 

P3 XZ - {YX + Y-^X)Z^ + XZ^ 

Pi = XZ-{1+XY)Z^ + YZ^ 

The chiralities {N, E, S, W} and the location (i, j) are simultaneously measurable, so the probability 
of a QRW started at ((0, 0), N) to be found at (r, s) at time t is the sum over 1 < fc < 4 of \Cr,s,tPk\'^ , 
the squared moduli of the probability amplitudes of going from ((0, 0), N) to ((r, s), ^) in time t for 
^e{N,E, S, W}. 



A number of different two-dimensional QRW's are analyzed in BBBPO^. In these examples 



the variety defined by the common denominator det(/ — M) of the entries of (/ — M) ^ turns out 
to be smooth, and amplitudes may be computed f rom Theo rem We now restrict our attention 



to the Hadamard QRW. This is not analyzed in |BBBP08[, because the denominator factors, one 
factor having a quadratic point at ±(1, 1, 1). 



The amoeba and the normal cone 



Denoting Hi = 1 — Z,H^i = 1 + Z, we have Fj — Pj /{QHiH^i), which is in the format of the 
quadratic point hypotheses with rj — 2. The origin is on the boundary of a component B of the 
complement of fiamoeba((5i?i^f-i) containing the negative z-axis. As before, the cone tano(i?) is 
the intersection of components of fiamoeba((3)^, fiamoeba(_ffi)^, and 6amoeba(iJ_i)'^. The latter 
two are just the halfspaces {{x,y,z) : z < 0}, which are equal and contain the component Bq := 

(a; + 2 + y^)/2, which we recognize from the Aztec Diamond example. 
Therefore, B ^ Bq and B* = {t^ > 2{r'^ + s^)} as in Section [41] 



Classification of critical points 



The intersection of Vq with Vh^ is the set 

^ . x + x-i + r + y-i 
Z ^ 
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The varieties Vhi and Vh_i do not intersect. We therefore stratify by 

Vi := {(1,1,1), (-1,-1,-1)} 
V2+ := VqCiVhA^i 

V3+ := Vh, \(ViUV2+) 
V3- := Vh_, \(ViUV2-) 
Vi Vq\(ViUV2±UV3±) 

Again, we are interested in the region of non-exponential decay, where x^ax = 0; checking where 
the strata intersect the unit torus, we find that the strata V2± do not intersect the unit torus. The 
strata V3± intersect the unit torus on the set {\X\ ~ \Y\ = Z = 1}, and the logarithimic gradient 
is always in the t direction. The factors Hi and H^i cause this direction to be obstructed and we 
are therefore not able to say anything about asymptotics in the (0,0, t) direction. This direction 
corresponds to the bright spot in the middle of the amplitude intensity plot in figure [TBI where there 
appears to be a bound state (probability amplitude for being precisely at the origin does not decay 
with time). 



Contributions from Vi occur for r interior to N*. The big-0 estimate, Lemma 1531 below, allows 
us to bound the magnitude of these contributions. We first compute the homogeneous degree of F 
at the point (1, 1, 1). The factor 1/Q has degree —2 here, the factor 1/Hi has degree —1, and the 
factor I/H2 has degree zero. The numerators Pj vanish to order two at (1, 1, 1) for all 1 < j < 4. 
We therefore have deg(F, (1, 1, 1)) = deg(/, (0, 0, 0)) = 2 — 2 — 1 = —2. Applying the lemma, we find 
that for the quadratic point w = 0, we have 6contrib(w) = 0{t^^). 



Finally, we compute the contribution from V4. It is shown in BraOTj ] that for each r interior 
to N*, there are precisely four smooth critical points on the unit torus, a conjugate pair and its 
negative: Z(f), Z(f), — Z(f), — Z(f). Denoting logZ(r) = z = iw, Theorem 13.41 tells us that 

6contrib(w) ~ C(f)|rp^ exp(— ir • w) 

where the magnitude of C(f) is proportional to the — 1/2-power of the complex curvature of log Vg 
at z. Adding this to the contribution from Z, namely 6contrib(?i;z/;), we obtain a quantity whose 
magnitude is 2cos0(r) times the magnitude of 6contrib(w), where 6 is the argument of 6contrib(w); 
note that differs from r • w by 7r/4 because the curvature is complex and its —1/2 power has 



argument — a7r/4, where a is the signature of q, which in our case is 1 — 2 = —1; see [BBBP08 , 
Section 2.3] for details on the phase of the curvature. Adding the contribution from the negatives 
of these two points kills the terms for which r + s + t is odd and doubles the even terms. This 
corresponds to periodicity of the walk. For fixed i, the phase term cos0(r) varies rapidly (with 
period of order 1). Ignoring the Moire pattern resulting from this term, the probabilities are of 
order and are spread over the disk + = which is the slice of the normal cone at the 
fixed value of t. 
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Figure 16: two-dimensional Hadamard QRW probability amplitudes at time 200 



4.4 Friedrichs-Lewy-Szego graph polynomials 

In the the study of a discretized time-dependent wave equation in two spatial dimensions, Friedrichs 
and Lewy required a nonnegativity result for the coefficients of whe re Q{X, Y, Z) :— {1~X){1 — 
Y) + {1-X){l-Z) + {1-Y){l-Z). To solve this problem, Szego [Szesd ] showed that the coefficients 
of Q^^ arc nonnegative for all /? > 1/2. Scott and Sokal [SS06| later observed that Q is a special 
case of a spanning tree polynomial of a graph. They proved a generalization of Szego's result to all 
series-parallel graphs. Their results are proved via the stronger property of complete monotonicity 
and are related to the half-plane property. In order to investigate whether these results might hold 
for the polynomials of a larger class of graphs, Scott and Sokal needed a means of checking the 
asymptotics of the coefficients: asymptotic nonnegativity is a necessary condition for term- by-term 
nonnegativity. 

The simplest nontrivial case in which asymptotics may be worked out is the one above. Szego's 
1933 proof of nonnegativity was, according to Scott and Sokal, "surprisingly indirect, exploiting 
Sonie-type integrals for products of Bessel functions." It is evident that asymptotics in this case 
may be derived directly from Theorem l3.7l We remark that the connection between these coefficients 
and harmonic analysis of symmetric cones is known to Scott and Sok al, who exploit the connection 



and cite several results on the subject from the sources FK9J; |lshOC | . 



We first check that is on the boundary of the component B of fiamoeba((5)° corresponding to 
the ordinary power series 1/Q = J2r c^rZ"". This follows if we show that Q{X, Y, Z) ^ for X, Y, Z 
in the open unit disk. To see this, let Di denote the open unit disk, let D2 denote the open disk 
{\z — 1| < 1}, and let D3 denote the halfspace {z : Re{z} > 1/2}. The set D3 is the image under 
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z 1/z of I?2 and D2 = 1 — Di. Therefore, Q has a zero on the open unit polydisk Df if and only 
if XY + YZ + ZX has a zero on Z?!; this is equivalent to 1/X + 1/Y + 1/Z having a zero on 
which is equivalent to X + Y + Z having a zero on D^. This is impossible because is contained 
in the open right half-plane. 

Composing with the exponential, then taking the leading homogeneous part, we obtain 

q = fihom((3 o exp, 0) ^ xy + xz + yz . 

We recognize this as half the quadratic factor in Section [42l Therefore q* is twice what is was there: 

q*(r,s,t) = 2{rs + rt + st) - (r^ + s'^ +t'^) . 

Let P{Z) = 1. Recalling that M is chosen so that the matrix for the quadratic form is (M^^)^ DM^^ , 

1 1 



2 2 

1 i 

2^2 

i i 
22" 



for q we obtain 



we see that the determinant of A/ is det(g) plugging in the matrix 
\M\ = 2. Thus, for /3 > 1/2, equation (|53)) gives 

as r varies over compact subsets of the cone 2(rs + + st) > + + i^. To check this estimate, 
let r = s = t = 50 and compute 

41-/37500/5-3/2 

Or ~ ^ ,^ w 0.000222832 . . . 

V^r(/3)r(/3~i/2) 

when (3 = 3/4. We then use Maple to crank out the true value of 059,50,50 which is rational mmrbcr 
near 0.000223464, for a relative error of around 1/400. 



4.5 Superballot numbers and multiset permutations 

Gessel |Ges92| defines the super ballot numbers by 

, , , (fc + 2r)!(2n + fc-l)! 



(fc- l)!r!n! (n + fc + r)! 

k f 2n + k\ 

These are a generalization of the ballot numbers ( ) (obtained by setting r = 0), 

2n + k \ n ) 

which are in turn a generalization of the Catalan numbers (set fc = 1). The Catalan number and 
the ballot numbers are integral and have combinatorial interpretations. Gessel shows that the super 
ballot numbers are integers as well and sets as a goal to find a combinatorial interpretation. 

After re-indexing via -B(a, 6, c) := 3(0, 5 — a — c, c) for 6 > a + c, one may extend this definition to 
all nonnegative (a, 6, c) and obtain the generating function F{X, Y, Z) — ^ B{a, b, c)X°-Y''Z'^ 
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from equation (|1.7p . The numbers B{a,b,c) satisfy the same recurrence as the nu mbers N(a, h, c) 



defined to have generating function 1/(1 — X — F — Z + AXYZ) and shown in Ask75l | to have 



nonnegative coefficients. There is no direct combinatorial interpretation know n for N(a, b, c) but 



it may be represented as a difference of cardinahties of multi-set permutations [AI76|. Gessel goes 
on to find several more identities involving these numbers and their generating functions, but no 
asymptotics are derived. 

The coefficients of F are of greater interest than the coefficients of G, but the fractional power on 
the non-quadratic term takes this problem beyond the main results of this paper. The deformations 
in Section refsec:homotopies still apply, but the further analysis in Section [^21 via Leray cycles does 
not work when this factor is algebraic rather than a simple pole. We therefore do not state detailed 
asymptotics for reserving this for future work. 

Let Q{X, Y, Z) -.= 2 — X — Y — Z + XYZ, so that 2/(3 is the ordinary power series generating 
function for 2~'^^''^^A^(a, 5, c). It is easy to check (e.g., via Grobner bases) that Q and its gradient 
vanish simultaneously exactly at the point 1. We have q :~ Q o exp = xy + yz + xz + 0{\{x, y, z)\^, 
whose homogeneous part (at 0) is given hy q — xy + yz + xz. Again, 

q*{r, s, t) = 2(rs + rt + st) - (r^ + s'^ + f) 

and \M\ = 2. 

There is a component B of fiamoeba((5)^ containing a translate of the negative orthant (cor- 
responding to the ordinary power series expansion); let us check that G dB. Proceeding as in 
Section it suffices to verify that Q has no zero in the open unit polydisk Di, which is equivalent 
to checking that Q{1 + Z) has no zero in — where D2 = {z : \z + 1\ < 1}. We have 

Q(l + Z) = XYZ + XY + XZ + YZ = XYZ ( 1 + — + — + - 

^\ ^ J -r -r -r \ X Y Z 

whence this is further equivalent to 1 + X + Y+Z having no zero in —D^, where —D^ is the halfplane 
{z : Re{z} < —1/2}. This is obvious, because the real part of X + Y + Z is bounded above by —3/2 
on 

We now apply CoroUarv 13.81 to obtain the asymptotics of 2~'^^^^'^N{a,b,c) inside the cone N*, 
these asymptotics being exponentially small outside N*. Letting P{Z) = 2, we plug P, \M\ and q* 
into p.7p to obtain the leading term asymptotics for the coefficients of 2/Q, which give 

N{a, b, c) - 2'^+^+'= — (2a6 + 2ac + 2bc - - b^ - c^Y^I'^ 
2tt 

uniformly on compact subcones of N*. For example, if a = 1, 6 = 20, c = 30, then the approximation 
yields N{a,b,c) w 2.595 x 10^^ while the actual value of iV(10, 20,30) to three decimal places is 
2.547 X 10^6. 
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5 Homotopy constructions 



Recall from the heuristic discussion following (jl.Sp that moving the chain of integration in (|1.2p to 
the torus ReLog~^(xinax) is not enough. Our goal in this section is to construct homotopies moving 
this chain of integration, within the domain of holomorphy of the integrand (but not within the 
domain of convergence of the Laurent series) to a different chain on which the maximum modulus 
of the inte grand is small except in a neighborhood of ficrit(r ). While the Morse theoretic methods 



GM88l | are in principle constructive, we follow ABG70|, taking advantage of hyperbolicity in 



of 

order to produce vector fields along which chains may be shifted. 



5.1 Vector fields 

Let Tr := {R/2'kIjY denote the d-dimensional flat torus. Given a Laurent polynomial F and a 
component B of the complement of 6amoeba(i^) , pick a unit vector f in the interior of the convex 
dual —B*. We suppose that r is a proper direction for B. Let 

U= \J (5.1) 

we6W(r) 

be the disjoint union of neighborhoods of each w G fiW(r), where 6W(r) arc the logarithmic critical 
sets from Definition 12.211 

Lemma 5.1 (Vector field away from the critical set). Let F he a Laurent polynomial, let B he a 
component of \ fiamoeha{f) , and let r he a unit vector in the interior of the convex dual ~B*. 
Suppose f • X is maximized at a unique x„iax in dB and define the neighhorhood U of flW{r) as 
in (|5.ip . Then there is a smooth vector-valued function rjjjo : Tr \ J7 ^ M'' such that: 

(i) iiviy) e K-^'-^(exp(xmax + «y)); 

(ii) f • Tju^ (y) = 1 for all y G Tr \ [/. 

Proof: First, for each y ^ [/, we will find a neighborhood Afy and a vector Vy such that ryj/c = Vy 
fulfills {i)-{ii) on My. 

Fix y . If /(xmax + iy) 7^ then choose a neighborhood My of y in such that for v S My, 
the quantity /(xmax + *v) does not vanish. Choose Vy with f • Vy = 1. 

Alternatively, suppose that /(x„iax + *y) = 0. By Proposition 12. 121 the homogeneous part, call 
it Ay, of the function v i— > /(x^ax + iy + v) is real and hyperbolic, and by Proposition 12. 81 there is 
a cone K of hyperbolicity containing tanx,„^3,(-B). Also by the first part of Proposition l2.221 there is 
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some Vy £ K with f -Vy = 1. By semi-continuity (part (i) of Corollarv l2.15p . Vy G K(cxp(xmax+*u)) 
for every u in some neighborhood My of y. 

The coUection {A/'w : w ^ U} covers V^] shrink it shghtly if necessary so that the closure of its 
union does not intersect 13W(r) . We may (shrinking some of the rj^ shghtly if necessary) choose a 
finite subcover {A/'w : w G 2} whose union has closure disjoint from fiW(r). Choose a partition of 
unity {ip^ : w G S} subordinate to the subcover. Define 

m-iy) ■= X! ^w(y)vw(y) ■ 

Now [i) is satisfied by convexity and (m) is satisfied by linearity. □ 

Corollary 5.2 (Vector field defined everywhere). Under the conditions of Lemma \5.1\ there is a 
smooth vector field rj satisfying 

7j{y) G K/'^(exp(x + zy)); (5.2) 
f • 77(y) = 1 onU". (5.3) 

Proof: Let ryw : U^r C'* be any map for which ?7w(y) G K^''^(exp(x + iyj). To see that we 
may choose such a map smoothly, note that the constant map ?7w(y) = v is such a map whenever 
V G tanx(i?). The reason for allowing a general function rj^ in place of a constant is that later we 
will use (|5.4p with functions rj-^ tailored to more specific needs. The collection {A/'w '■ w G S}U {[/w : 
w G fiW(f)} covers Tr. Choose a partition of unity {'(/'w} subordinate to this and define 

weH wGfiW(f) 

This proves the corollary. □ 

We remark for later use that if f is replaced by a non-unit vector r, then applying the above 
constructions to f replaces (15. 3p by r • ry(y) = |r| on If^. Next, we give a projective version of the 
above construction. We say that a 1-homogeneous function is smooth if it is smooth away from 
the origin. 

Lemma 5.3 (projective vector field). Let A be a real homogeneous polynomial in d variables of 
degree m > 1 and let B be a cone of hyperbolicity for A whose dual —B* has nonempty interior. For 
each y G M.'^, recall the cone K"^'^(y) defined in Proposition \2. 8\ Let r be a non- obstructed vector 
in the interior of —B*. Then there is a 1-homogeneous, smooth vector field rj on M'^ such that for 
all y G K'' and all r' in a neighborhood of r, 

ft) rj{y) G K^^^(y); 
fti) r'-rj{y) > |r'||y|. 
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Proof: This is a homogeneous version of the proof of Lemma [5. II Assume first that |y| = 1. Wc 
define 77 locally and then piece these together via a partition of unity. When A{y) ^ we can find 
neighborhoods J\fy of y and Af* of r such that A vanishes nowhere on Afy and there is a v for which 
r' • V > |r'| on Af*. By the trivial part of the definition, Vy € K'^'^(y). 

When A{y) = 0, because r is non-obstructed, there is a vector Vy e K"^'^(y) with r • Vy > 0. 
By semi-continuity (part (m) of Corollarv l2.15p . Vy G K"^'^(u) for every u in some neighborhood 
Afy of y. By continuity, r' • Vy > for every r' in some neighborhood JV* of r. We may then replace 
Vy by some positive multiple so that r' • Vy > |r'| for r G J\f*. 

To define the 1-homogeneous function ry, it suffices to define it on the set 5*1 of vectors y of 
norm 1. Cover Si by finitely many neighborhoods {A/'w : w G S} and use a partition of unity 
subordinate to the cover to define 77 via (|5.4p on S'l. Extending 1-homogeneously via 77(Ay) := Ary(y) 
finishes the construction. □ 



5.2 Homotopies 

Given any triangulation of the torus Tr by oriented d-simplices, one may regard any continuous 
function </> : Tr ^ as a d-chain on f2 (composing with the singular simplices of the triangulation) . 
The integral of d-forms on such a chain is independent of the triangulation, so we may refer to cf) 
itself as a d-chain. Let 77 be any continuous vector field on Tr and fix any e > 0. Define the homotopy 
$ = : Tr X [0,1] ^ by 

My) := 7y + X + e [(1 - t)u + tr]{y)] . (5.5) 

We specialize now to 77 given by Corollarv l5.2l Setting t = gives a cycle whose range is the torus 
Tu := — eu -I- iTr. Setting t = 1 gives another cycle, which we call C,,, shown to be homotopic to T 
in via the homotopy {$t:0<i<l}. 

Theorem 5.4 (The homotopy defined by 77 avoids V/). Let -q satisfy (|5.2p - (l5.3p and define $t 
hy (|5.5p . Then for e > sufficiently small and all <t < 1, /($t(y)) 7^ 0. 

Proof: This follows from Corollarv l2.16l □ 

Let {C/w : w G fiW(r)} be disjoint open neighborhoods of the points of 6W(r) as before, and let 
Cw denote the restriction of to U-w- Let Cjj^ denote the restrition of to C/^. The chain C,; is 
representable as a sum 

Cjj = Cu" + ^ Cw • 

wGfiW(r) 

An immediate consequence of the previous constructions is: 
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Corollary 5.5 (localization of the Caiichy integral). The chain is homotopic in (C/(27rZ))'^ \ V/ 
to a sum of chains 

where Cu-' and each dC^ are supported on the set ofz such that r-Re {z} — r-x > e|r|. Consequently, 
the integral (|2.ip decomposes as 



STTiy J^+^T. /(z) 



wej8FK(r) 



1 



e-''-^-—dz 



where R — 0(e~''""'), anrf Theorem \3.3\ follows. □ 



While the above general construction, e.g., with r]^ = u, suffices to localize the Cauchy integral, 
explicit computations for a cone and a plane, done in Section [6. 6[ will require specific choices of r]^. 
Say that tJw is projective if /^^(w + •) is homogeneous of degree 1 and smooth away from 0. The 
first of these two results follows immediately from Lemma 15.31 

Theorem 5.6. Let A be a hyperbolic homogeneous polynomial with cone of hyperbolicity B. Let r 
in the interior of —B* be non- obstructed and let rj be the projective vector field of Lemma \5.S\ . Let 
{<i>t} be the homotopy on defined by ((53)1 with x = 0. Then A{<^t{y)) ^ for all < t < 1 
except when t — 1 and y — 0, andr-$i(y) > c|y|. Consequently, for u £ tanx(-B), the chain u + iM.'^ 
is homotopic in the complement of Va to the projective chain Cq :~ <i>i[M'^] on which r • y grows 
linearly in |y|. This construction is uniform as r varies over some neighborhood M . □ 

Definition 5.7. Let Cq denote the chain parametrized by Tr defined by 

y^^y + e[(l - t(l -{6- |y|)+)))u + t(l - i(l - (S - |y|) + ))r;(y)] . 

This is obtained by replacingt is the definition o/Cq byt{l — {S—\y\)~^); the definition does not change 
the homotopy outside the ball of radius 6, but inside this ball the homotopy stops early, stopping at 
time 1 — 6 at the origin and interpolating linearly in \y\. 

Theorem 5.8 (locally projective homotopy). Let f be a Laurent polynomial and let B be a compo- 
nent ofW^\fiamoeba{f). Suppose r is proper with dual point Xmax, that flW{r) is finite, and that r is 
non-obstructed. For each w G flW(r), let f^ ~ phom{f,'w) and let r/w be the vector field constructed 
in Lemma [573\ with f^ in place of f. Let {C/w : w e pW{r)} and U be defined as in ()5.ip . define r] 
by ([Q]) . and define {^Y'} by (|531l . 

Then, if e > and the neighborhoods {C/w} are taken to be sufficiently small, the homotopy 
will avoid V/. Also, the inequality 

r' • *i'i(y) - r' • X > (5 min |y - w| (5.6) 
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will be satisfied for some S > 0, for every y S Tr and every r' in some neighborhood of v. 



Proof: By consruction, Vw(y) G K-''"^-^(y). By semi-continuity (part (i) of Corollary 12. ISp . also 
Viu(y) G K^'^(y) for y in some neighborhood U!^ C U^. Replacing the sets ?7w by sufficiently small 
sets, we therefore have Vw(y) G K-^'^(y) for all y G U^. The argument is now the same as the proof 
of Theorem [El □ 



5.3 Consequences and an example 



As outlined in Section [TT2l the deformations constructed in Theorems l5.4l and l5.8l allows us to localize 
and then compute the Cauchy integral. These computations are carried out in the next section using 
Fourier ap paratu s. We record here a preliminary estimate that is useful in a more general context 



(see, e.g., [Bra07|). If 

k 

F (5.7) 

is the product of d-variatc polynomials to arbitrary real powers and Z is any complex vector, the 
homogeneous degree deg(F, Z) of F at Z is defined by 

k 

deg(F,Z) :=^s,dcg(Q„Z) 

(it is easy to check that this is independent of the representation of F as such a product). 

Lemma 5.9 (big-0 estimate). Let F = 11^=1 Q^j '^"■^ ^ denote the product of all the Qj for 
which j is not a positive integer, so that Vq is the singular locus of F. Let f := F o exp and let 
B be the component of \ fiamoeba{G) containing a translate of the negative orthant (whence the 
corresponding Laurent series is the ordinary power series expansion). Fix a proper, non-obstructed 
direction r in the interior of —B*. and let be the locally projective homotopy constructed in 
Theorem 1 5. 8\ via the neighborhoods and vector fields rj^ for w G pW(r). Then 

(i) Lf(j){z) = 0{\zf) and (3 + d>0 then 

Iz"] J exp(-r-z)(/)(z)rfz = 0(|rr'*-^) . (5.8) 

(ii) Consequently, i/deg(/, w) + d > 0, then for any bounded funtion ip, the following integral is 
absolutely convergent and 



j e-'-Xz)/(z) dz = (^|r|-('*+d<=g(/.w))^ 
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(in) It follows further that the Taylor coefficients of F satisfy = O {\z\ '^\r\ ") where 

a — d + max deg(/, w) . 



Proof: The chain Cw is a cone avoiding the singular locus of the homogeneous function /. Let 
S denote the section {z : Rc{f • z} = 1} of this cone. We have Cw — [0, oo) x S. Write /(z) = 
\z\'^°^^^''^'> Fq{z/\z\) for some smooth function Fq on 5 and decompose dz — f^^^ dt A dS for some 
form dS on S. Let M :— Jg \Fq{u)\ dS{u.) and M' :— sup|-0|- Integrating first over S then over 
[0,oo) gives 



z / e 



V(z)/(z)dz 



^dog(/,w)^d-l^^g-|r|t 



i^{u)Fo{u)dS 



< 



g-|r|t^d+dog(/,w)-l MM' dt 



which is absolutely convergent and 0(|r 



-(d+dcg(/,w)) 



) as desired. 



We have seen that the chain (x + u) + zTr is homotopic to the sum Cw + X]we6W(r) 
domain of holomorphy of /. Consequently, the Cauchy integral 



/ e-"-"/(z)dz 

Jx+u+iTr 



is equal to the sum of the integrals over Cw and Cw- The last result now follows. 



□ 



We close the section with an example of the conclusion of Lemma 15.31 and Theorem 15.81 in the 
case of the product of a quadratic cone Q and a linear function H . We give an explicit construction 
a vector field rj satisfying the conclusion of the lemma; this explicit construction will be useful in 
computing an integral in Section 16.61 

We consider the case of F = 1/{QH), where q :— Q o exp and h :— H o exp are respectively 
quadratic and smooth at the origin: 



q 

h 



q + Ri 
h + R2 



with q homogeneous quadratic, h linear, i?i(y) — 0(|yp), and i?2(y) = 0(|y|^). The signature of g 
is assumed to be (1, 2) whence the zero set of q in real space is a cone over a circle. Suppose that 
the zero sets of q and h intersect transversally in real space. In other words, the plane {h = 0} 
intersects the cone {q = 0} in two lines; in projective space, the line {h = 0} intersects the circle 
{q = 0} in two points. 

The construction of rj depends on the choice of the cone of hyperbolicity B of qh; in the applica- 
tions below, this is the component of fiamoeba(Qi/)'^ containing the negative half of the z-axis. Let 
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us assume in this example that A := qh is hyperboUc with respect to —63. We also assume without 
loss of generality that h{—e^) > 0. We have seen in several of the examples that B :— Bi ni?2 where 
Bi is a halfspace dual to h and B2 is the projective ellipse defined by q. The dual cone is the cone 
over the teardrop pictured in figure 13.31 here the conic is the dual to Bi , and the vertex is the line 
dual to the hyperplane B2', as usual, the dual to the intersection of the convex sets is the convex 
hull of their respective duals. 

We will construct the section rj guaranteed by Lemma 15.31 for which r • ri{y) > 0, along with a 
null section 77 satisfying r ■ f] = that is needed in Section 15.31 

Fix r G N*. There are two nearly identical cases, depending on whether or not r G -62- Assume 
first that r G -B|. In figure [T71 below, height is the linear functional defined by — r, so the plane 
Xr := {r • X = 0} is drawn as as horizontal, with r • x increasing as one goes downward. This plane 
intersects the real cone {q — 0} only at the origin, hence the cone has a positive (lower) and a 
negative (upper) half; we have assumed —63 is in the upper half. The construction of 77 is automatic 
if we mandate that fi{y) = — esjyl + Ay for some A. In order to obtain r • ^(y) — 0, we need to take 

A = |y|^. (5.9) 

Wherever r • y ^ 0, this is clearly smooth and 1-homogeneous. Setting aside the points where 
r • y = 0, at every other point y where q ovh vanishes but not both, the cone K"*'^(y) is a halfspace 
bounded by the tangent plane at y to {A ~ 0}. This halfspace contains the vector y, making it 
obvious that of the two halfspaces bounded by this plane, f}{y) is in the one containing —63, thus 
is in K"^'-^(y). When y is in the intersection q = h ~ 0, again ?7(y) G K.^'^{y), because this cone 
is the intersection of two halfspaces, each of which we have seen to contain ry(y). Finally, to deal 
with the points where r • y = 0, note that q is nonvanishing (outside of the origin) on this set. In 
a neighborhood of the point ft, = 0, we use a smooth bump function ip^ : M. ^ [0, 1] that is one on 
[— e, e] and zero outside [— 2e, 2e]. Letting x be any vector with x • r = and ft.(x) > 0, we take 

viy) ■■= \y\M^ ■ y)x + (i - A) {-\y\e3 + Ay) 

where A is defined by (|5.9p . This completes the construction of 77. When ^(y) ^ 0, the cone K^'^(y) 
is all of K'^, so verification that ri{y) G K"^'-^(y) is trivial; we conclude that 77 is a I-homogeneous 
section of K'^'^(-) with r • 7/ = 0. 

What we have accomplished is to find a single formula for rj that works for all strata of {A — 0}, 
resorting to partitions of unity, only in one place, away from {q = 0}; this will be useful in the 
sequel. For a pictorial description of the construction we have just completed, see figure [T71 In the 
upper half of the cone, rj points inside, as does —63. In the lower half (not shown), both fj and —63 
point outside. To obtain the vector field of Lemma l5.31 we set 77(y) = fi{y) +e|y|e3 for a sufficiently 
small e. 
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Figure 17: the vector field 77 



When r is outside the circle, but inside the teardrop, the plane orthogonal to r does intersect 
the real cone q = 0. In projective space, the analogue of figure [T7] is figure [TS] This case is in 
some ways simpler because one may choose a vector v in the cone B2 for which r • v — 0. Because 
B2 C K'^'^(y) for every y where q vanishes, setting ?/(y) = v works everywhere except where h 
vanishes and K"^'^(y) may be smaller. In a neighborhood of this projective line, we may instead 
take fj{y) = — y — ce^ for some c > 0. Piecing these together, projectively, via a parition of unity, 
finishes the construction. 



1 = 




li = >■ = (' 

Figure 18: when r is outside B2 

Finally, we note that when r is on the dashed boundary in figure 13.31 it is obstructed and the 
above construction does not work. 
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6 Evaluation of integrals 



6.1 Reduction to integrals of meromorphic forms over projective cycles 

In this section we prove Theorems 13.71 and 13.91 The first step is to show that the hnearization in 
Lemma 13.21 may be integrated term by term, and also that the chain of integration may be replaced 
by the cone Cg. With Cw as in Theorem 15. 8[ we recall from p.3p the quantity 

ficontrib(w) := ( — \ f e"'"-^ ^ dz . 

K^^iJ Jc„ q{^yU,=i^i 

Lemma 6.1. Let F satisfy the quadratic point hypotheses. Let Cq be as given in Definition \5. 7[ Let 
c(m, /, n) be the coefficients of the expansion given in Lemma \S.'2\ for the function f := F o exp at 
the point x^ax + Then for every N > 1 and sufficiently small i5 > 0, 



^flcontrtbi^) - f ^cim,£,n)y'^qiy)-^-'f[h,{yr^--dy (6.1) 



(2TTi) 

^ ' |m|-/i<'-fcn<JV^O j=l 



In the case k — 0, this reduces to 
Z 



(27r«)° 



^-^flcontrtb{yv) - V / c(m, ^)y'"q(y)-^-' dy + O (|r| 



2s-d-N\ 



Proof: Let C be the common restriction of Cw and logz + Cg to a neighborhood where they agree. 
By (|5.6p . if 6 is sufficiently small, we may truncate Cw to C losing only an exponentially small 
quantity: 



{2mf ficontrib(w) 



exp(-r • y)/(y) dy 



exp(-r • z)/(y) dy + Ooxp(exp(-r • z)) 



-logz 



exp(-r • y)/(y) dy 



-'exp 



(1) 



Use Lemma [3.21 to write / as a finite series plus a remainder term. Integrate term by term over C 
and apply the big-0 estimate (|5.8[) to the remainder term to get the conclusion of the lemma with 
Cq replaced by C — logz. Replacing C — logz by Cq affects the integral by Oexp(l), estabhshing the 
lemma. □ 
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6.2 Generalized functions 



The integrands in (|6.ip are hom ogeneous and the chains of integration projective. Many such 



integrals are evaluated in ABG70l | but there is a hitch: the integral is evaluated not over Cq but over 
iW^. The latter integral is in general not convergent over iM."^ for two reasons. First, integrability will 
fail near the zeros of q whenever s is large. Secondly, because | exp(r • x)| = 1 on iR'' (as opposed to 
the exponential decay on Cw), integrability at infinity will fail whenever |m| > 2s— d. These problems 
are solved respectively by moving the contour and by inserting compactly supported functions inside 
the integral. The apparatus to do this is the theory of generalized functions (distributions) and their 



Fourier transforms, developed in GS64| and elsewhere. We summarize the facts needed from this 
literature. 

We work with two linear spaces that are dual to each other. We call these M"^ and M"**. We fix 
bases dual to each other so that for r e W^* and x e R'', we have {r,x) = X]j=i '^j^i- Hereafter, 
we denote the pairing by r • x. While all of the ensuing constructions could be defined on either 
space, our purposes require slightly different constructions the two spaces and we reduce confusion 
by developing these asymmetrically. 

Let Co(R'^*) denote the space of s mooth complex valued functions on R''* with compact support. 
These are called test functions in |gS64| and the closed support of a test function g is denote 
by 6supp(g). Topologize test functions by convergence of all derivatives; this may be metrized, for 
example, by 

n |k|=n ^ 



where (/)(x) — x/{x+l). The space G* of generalized functions (sometimes called distributions) 

is defined to be the dual of Co(R'^*), namely the space of continuous linear functions on Co(R''*). 
Let filoc-int be the space of locally integrable functions on M.'^*, that is, functions g such that 
g G L^{Bj\j) for the ball Bpj of every radius N in R'**. There is a natural embedding of filoc-int 
into Q* mapping the function / to the linear map g ^ j f{r)g{r) dr. We denote by t/ the image 
of / under this identification. Generalized functions in the image of this identification are called 
standard functions, but there are many nonstandard functions. One example is the function Sr 
defined by Sr{g) = g{r). 

Sometimes a function is not standard but agrees with a standard function on some region. Let 
V be an open set in (R'')* and suppose that for any g whose closed support is contained in V, the 
value of the generalized function L is given by f{v)g{r) dr for some function / G L^{T>). We then 
say that L is partially identified with / on T). 

Differentiation may be defined on Q* by 

^ C:^g^-C (^g\ . (6.2) 
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This commutes with the identification map: integrating by parts, 

= J A/(r)5(r)dr (6.3) 

which is evidently the embedded image of df /drj applied to g. An example of this is the generalized 
function {d/dri)6r which maps g to {dg/drj){r). A famous result (not needed here) is that every 
generalized function is of this form: given C E G*, there is a continuous / S filoc-int and a k for 
which C — d^f /dr^. Restricting the integral to V, we see that differentiation also commutes with 
partial identification. 

On M'' we define a slightly different space of test functions. Denote by CRD(IIi'') the space of 
rapidly decaying smooth functions, meaning that they are 0(|x|^^) at infinity for every > 0. 
Again, topologize by convergence of all derivatives. Let Q denote the dual of Crd(K''). This space 
of generalized functions is slightly smaller than the space Q*. Let 6poly-bd denote the space of 
functions / on satisfying |/(x)| < C(l + |x|)^ for some C,N > 0. Then the space fipoly-bd 
embeds in Q; again, we denote the image of / under this identification by if. 



6.3 Inverse Fourier transforms 

We now define Fourier transforms and their inverses. Fourier transforms will be defined for functions 
on the dual space, while inverse Fourier transforms will be defined for functions on ordinary space. 
Fourier transforms will be defined only for nice functions, while inverse Fourier transforms will be 
defined for generalized functions. 

For g e Co(M''*), define the Fourier transform g by 

5(x) := / 5(r)exp(-ir-x)dr. (6.4) 

Observe that g G Crd(R'') (this is the Riemann-Lebesgue Lemma). In fact, we may extend g to 
a f unction on all of C''. This is a holomorphic function and for every integer Af > it is shown 
in |ABG7Cll . (2.3)] to satisfy an estimate 

|g(x + zy)|<C(Ar)(l + |x + zy|)-^exp[ sup r-y). (6.5) 

\refSsupp (g) / 

Let >C be a generalized function in Q. We define the inverse Fourier transform J-^^{C) by 

:F~\C){g):={2nr''C{g). (6.6) 
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This is well defined because we have just seen that g G Crd(K'^) and it is easy to see that it is 
continuous and therefore an element of Q* . Suppose that / G L^iW^). Then 

T-\if){g) = {2ir)-\f{g) 

= (2^)-'^/ /Wf/ 5(r)exp(ir.x)dr') dx. 

Since |/| and \g\ are integrable, we may switch the order of integration to see that J-^^{if ) is the 
generalized function identified with the actual function 

i2nr'^ J /(x)exp(ir •x)dx. 

Although we do not need it here, we remark that inverts the Fourier transform: let g e CoiMf^*) 
so that g S L^{M.'^); then the above computation shows that T~^{g) is equal to the standard function 
(27r)~'* J g(x) exp(ir -x) dx, which is equal to g by the usual theorem on inverting Fourier transforms. 



Boundaries of holomorphic functions 

To evaluate the integrals arising in this paper, we must examine generalized functions arising as 
limits of holomorphic functions. Let / be holomorphic in a domain M.'^ + iA where zero is contained 
in the boundary of A. Suppose that / satisfies an estimate 

|/(x + zy)|<C|yr~(l + |x|)^ 

for some iV > 0. The estimate (j6.5p shows that the integral 

f{x)g{x.)dx 

hi?) 

exists and is independent of G A. The same is true for J^^_^.^ f{x.)h{x.) dx as long as the esti- 
mate (16.511 is satisfied with h in place of g. In particular, this defines a generalized function in Q 



(see jABGTd . (1.3)-(1.5)] and following). We denote by la/ the generalized function this defines; 
if / has a limit in as 77 —> in A then oaJ is just this standard function. An example is the 
function /(x) — A(x)^'' for some homogeneous polynomial, A. If A is a cone of hyperbolicity for 
A, then / is holomorphic on R'' + lA and blows up no worse than a power of the magnitude of the 
imaginary part of the argument. When s is sufficiently large, this is not a standard function. 

Two classical and useful results generalize the analogous well known results for ordinary Fourier 
transforms. 

Proposition 6.2. Let f be a function satisfying /(x + iy) = 0(|y|^) for some N, as above. Let 
X™ be any monomial and let L be any linear transformation. Then the inverse Fourier transforms 
ofx"^f and f o L^^ are given respectively by 

am f-l( f] 

T-\x^f){v) ^ (6.7) 
T-\foL-')[r) = \L\T-\f){L*v). (6.8) 
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Proof: Pick any h e Co(]R'^*)((R"')*). Integrals in the following calculation will be over R'' + in 
the x-domain and over (M'*)* in the r-domain. Using the definition of Fourier transform in the first 
line, calculus in the second, integration by parts in the third line, Fubini's Thereom in the fourth, 
and integration by parts once more, we see that 



f x™/(x)/i(x) = f f x™/(x)/i(r)e*'"-^drdx 



^/»(^«^) h{r)drd^ 
//(r)(-4)'"Mr)dr. 



The left-hand side of this is {2ttYT ^(x™/)(/i) while the right-hand side is by definition (see (|6.2p ) 
equal to (27r)''ii'"[(a/5r)'"J^-i(/)](/i), thus verifying (l6J)l . 



The second assertion of the theorem is directly verified. Making the coordinate change x — Lx. 
and using r • Lx = {L*r) ■ x recovers (|6.8p . □ 

Suppose the function E on R''* is not locally integrable. Then lE is not a well defined generalized 
function. Nevertheless, lE is defined as a partial function: if 17 is a function supported on a set where 
E is locally integrable then i.E{g) = J E{r)g{r) dr is perfectly well defined. We wish to conclude 
that an acutal Fourier integ ral such as occurs in (|6.1|) of Lemma lOI is equal to the locally integrable 
function of r computed by ABG7Clj | as the generalized Fourier transform of the integrand in ()6.ip . 
We therefore require the following lemma. 

Lemma 6.3. Let F satisfy the quadratic point hypotheses and let Cq be as in Definition \5. 7| Let 

X™ 

m(x) :— ; — 

be any of the terms in the series expansion f at x + iw as in Lemma \3.Si Suppose that the inverse 
Fourier transform !F^^{u), defined relative to the domain — m -I- , is given by a partial function 
lT~^{u) on the set of non- obstructed dual vectors in the dual cone, N* to tanx(B). Let 

^(r) := (27r)-'' / e~''-yu{y)dy 

E{v) := tF-\u){-v). 

Then 

il^{v) = E{y) (6.9) 
for any non-obstructed r in the dual cone N*. 
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Proof: This is a matter of moving — iCp to — m + M** while introducing appropriate smoothing 
functions to maintain integrabiUty. We fix a neighborhood A/" of r of non-obstructed dual vectors in 
N*, as in the conclusion of Theorem 15.61 whose closure is in the interior of —B*. We then see that 







(6.10) 



exponentially fast in |y| as y — > cxd in — iCgi uniformly as r varies over Af. It follows that if 
9 ■■ (K'')* - 



is smooth and supported on some compact subset of M, then 
|^(y)| <c||5||exp(-c'|y|) 



(6.11) 



as y varies over Cg, where \\g\\ := J \g\ and c and c' are positive constants not depending on g. 
Now fix r G A/" and let gn be a sequence of smooth functions supported on Af and converging to 6r- 
Note that the estimate (|6.1ip holds for the function g ~ a,s well as for all g„, where in this case 
g{x) = 

To establish ()6.9|) . fix an e > 0. Nonvanishing of q and each hj on —iCq, together with (16. lip , 
implies that we may pick a compact set K such that 



\9n{y)u(y)\ dy < 



(6.12) 



for all n, and also for in place of g„. The sequence g„ converges to exp(ir • y) uniformly on K, 
hence we may choose A^o large enough so that for n> Nq, 



"(y) |exp(ir • y) - g„(y)| dy 



K 



< 



Increasing Nq if necessary, we may also ensure that 

E{r)-J E{v')g^{v')dv' 
for all n > Nq. We may now conclude that n > Nq implies 



< 



^(r) - {27:y 



9n{y)u{y)dy 



3 



(6.13) 



(6.14) 



(6.15) 



Indeed, the two terms we have subtracted are integrals over —iCq of two integrands; denoting the 
integrands by f3 and /3' , we break — iCq into K U K'^ and use the triangle inequality, viz.. 



< 



K 



and {2n)-'^ < 1 to obtain ((6T5)) . 



The homotopy —i^t in Theorem 1 5 . 61 moves —iCq to —iu + W^ while avoiding the singularities of u 
(by homogeneity, term is singular at y if and only if it is singular at iy) . We have seen in l|6.5p that 
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gn is rapidly decreasing on the image of this homotopy. Truncating at the boundary of a large ball 
and sending this to infinity shows that the integral in (|6.15p is unaffected by appliying the homotopy. 
Hence, 

(2^)-^ / .gn(y)w(y) dy - {2tt)-^ [ 9n{y)<y) dy = [ i?(r').9n(r') dr' 

by the definition of the inverse Fourier transform. This identity allows us to apply the triangle 
inequality to (|6.14p and (|6.15p . yielding 

\^p{r) - E{r)\ <e. 

Since e > was arbitrary, this proves the lemma. □ 



6.4 Proof of Theorem [37fl 



We begin with a result from ABG70|, evaluating the Fourier transform of 5 " where S is the 



standard symplectic quadratic 



For this special case, we let A := {x : a;i < 



and S'(x) > 0} be the cone of hyperbolicity containing the negative xi axis, we choose an element 
1] — —e\ of A, and we let N* := {r : ri > and — ^ 



fc=2 ' fc 



> 0} be the dual cone to A. In 
the case where s is not an integer, we also need notation to specify what is meant by S'(x)^''. To 
specify a branch of this is the same as specifying a branch of the argument function Arg5(x). On 
any simply connected domain where S" 7^ 0, this may be accomplished by specifiying ArgS'(x) at 
any point in the domain. Therefore, we write S'(x)~^|Arg(s(j)))=e to denote such a specification. 



Theorem 6.4 ( |ABG70l . Equation (4.20)]). Let S'(x) 



\ — I 2 ^ ■ ■ ■ ^ d- Then the inverse Fourier transform of S' 
s 0, d/2 — 1, it is given by 



c\ ~ x\ ~ . . . ~ x'ji, so that S*{r) = 
exists in a generalized sense and, if 



, s*{ry~w^'^ 

22.-i7r(<i-2)/2r(s)r(s + 1 - {d/2)) ' 



(6.16) 



To be precise, let A be the component of the real cone {S* > 0} that contains the negative xi axis 
and let rj G A, for example, rj = ei. Then if g is supported on a compact subset o/N*, 

(2^)-' / ^"^(x)g(x) f S*ir)gir)dr (6.17) 

where C = e"" / [2^''-^TT^'^-^y^T{s)T{s + 1 - d/2)] . □ 



Proof: This result is taken from ABG7Cll | but with definitions spread across several sections. So 
as to make th e citation checkable (especial ly in light of some minor errors), we reference a number 
of passages of ABG70|. Equation ABG7Cll . 4.13] defines a Fourier transform of S*"" (their notation 
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to 



Rie49|: 



for ig is g ). Then in ABG70l . 4.20] they give the following formula for this quantity, attributed 



S*{r) 



2~d/2 



T^d/2-i22s-iY(^s)T{s + 2- d/2) ■ 
The argume nt of the F function the second time is wrong: it should be s 



with 



(6.18) 

1 — d/2, agreeing 

Rie49|. They are also missing a factor of e*^*. To see that this factor should be present, note 
that their specification of the branch of is given at the top of page 146: they specify this over 
the simply connected set i-q + W'' by specifying that AYgS{iri) = tt + Arg {S{rj)). Taking 77 = ei, we 
see that Arg(S'(i7y)) must be an odd multiple of tt. For such a specification, the Fourier transform 
cannot be real for small real values of s. Indeed, taking s to be small and positive, and noting that 
S is never positive real on i'q + W'-, we see that all arguments of 5'(x) lie between 2nn and 2{n + l)7r, 
hence all arguments of S'(x)~'* lie between 2nsTi and (2n + 2)s7r. Let x = 177 + y. Then switching 
y and — y conjugates e^^'^ while reflecting 5'(x)~'' about the line Arg — {2n + l)s7r. Therefore, the 
integrand of the Fourier transform 



[ ue'^-'^Si^)-' dx 



is also reflected about the line Arg = {2n + l)s7r implying that the integral must therefore lie on 
the line of reflection. For small values of s this is not real, demonstrating the need for a correction. 
The corrected formula (|6.18p implies (|6.16p . □ 

Corollary 6.5 (Fourier transform of a cone). For any real quadratic A having signature (1, d — 1), 
any monomial x"", and any s ^ 0, d/2 — 1, the inverse Fourier transform o/x'"^"'* is given by 

\M\{d/dr)"'A*{ry~W2) 



22«-i7r(''-2)/2r(s)r(s + 1 - (d/2)) 
where M is any real linear transformation such that A = S o M^^. 



(6.19) 



Proof: Pick a linear transformation M such that A = S'o M~^. Recall that A*{r) = S*{L*r). Use 
the second part of Proposition 16. 21 and then the first to obtain (|6.19p . □ 

Proof of Theorem 13.71 We are required to prove the given asymptotic expansion of 6contrib(w). 
We have assumed no linear factors, so the second, simpler formula from Lemma l6.ll applies, giving 

ficontrib(w) = f'— Z-" f e-^ y ^ c(m, n)y™q(y)-''-" dy + O (|r|-(^+''-2s)^ 



as long as > 2s — d. 

Now integrate term by term. We see that 



ficontrib(w) = Z 



0(|r| 



-d+2s-N\ 



E 

|m|-2n<JV 



1-) I e-'-yc(m,n)y'"g(y)-^-" 



dy 



(6.20) 
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The specification of — s power for this generating function is that the argument of q(u) is zero. To 
turn this into a Fourier transform, the change of variables y = iy' is needed. Under this change of 
variables, dy = i'^ dy' , and the summand in (|6.20p becomes 



c{m,n){27T)-^ [ e'^^-^' iiyTii^y'r'-'' dy' , 



Now the argument of q{iy') is still continued from initial data Arg(g(u)) = 0, and this argument 
may also be written as in plus the argument of q{y') continued from initial data Arg m)) = — tt 
as y' varies over — iCq. The summand in (|6.20p now becomes 

c(m,n)(27r)-'^*l-le-"^ / e-^-''' lyT^V)"'"" ^y' ■ 

Everything is now lined up. Lcmma [6.3l shows that this integral, '(/'(r), is equal to the partial function 
E defined by the Fourier transform of the integrand relative to the domain — zu + K''. Corollarv l6.5l 
computes this partial function (recalling that the inverse Fourier transform builds in the factor 
(27r)^'*) and yields the summand 

1^1 |M|(9/9r)"'r(-r)-W^) 
cim, n}y i) 22«-i7r('i-2)/2r(s)r(s + 1 - (d/2)) ' 

Multiplying by the Z in front of the right-hand side of (|6.20|) establishes desired conclusion p.4p . 
still under the assumption N > 2s — d, which was used to bound the remainder term. 

Finally, if iV < 2s — rf then the foregoing argument may be applied with N replaced by the 
least integer N' greater than 2s — d. Each term in the sum with |m| — 2n > N is 0(|r|^'^+^'*^^) 
by (|2.12p : the remainder term satifies this bound as well because it is 0(|r|~''+'^*^^ ) with N' > N. 
The theorem is therefore proved for every N . □ 



6.5 Extra linear factors give rise to integral operators 

Let F satisfy the quadratic point hypotheses and let r be a non-obstructed vector in the dual cone 
to N* := tanx„,^(i3). Equation (16. 7|) of Proposition [672] has a moral inverse: just as multiplication 
by X turns into differentiation in the r-domain, division by a linear function in x should turn into 
integration in the r-domain. This subsection proves a theorem along these lines. However, because 
anti-differentiation is not well defined, the resulting formula (|6.2ip fails to specify which iterated 
anti-derivative will result. We show that the correct choice can be determined under the additional 
assumption 2s > 1. This is not, however, the case with the bulk of our examples, whence our 
alternative analysis in Section 

Let L be any linear function, with coefhcients i(x) — aixi + • • • + a^Xd- We may view L as a 

d Q 

vector in (M'')*. The notation d/dL will be used to denote the differential operator ^^Oj^ — on 
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9" / 5 \ 

r-space. Also, 'qJ~^ ■— Y\. ( qj^ ) denotes the corresponding sum of monomial operators (i9/9r)". 

Proposition 6.6. Let F — P/iQ'^ Yl'j^i H"' ) with rij positive integers, q quadratic and each hj 
linear. Let n — (rii, . . . , n^.) be the multiexponent of the functions Hi, . . . , in the denominator of 
F and denote p := n + nl = (rii + n, . . . ,nk + n). Then 

am 

ficontrtbiw) = ^"'■(-l)'P'+''"'c(m,^,n)— u,,„(r), (6.21) 

\m\-2t-kn<N 

where ui is an iterated anti- derivative of q^^^^ satisfying 



dhP 



ue,n = T-\q-'-') (6.22) 



gi7r(s+f)^|m| 



\M\{d/dr)'^A*{ry+^-W^'> 



22{s+e)-i^{d-2)/2Y(^s + £)T{s + £ + 1 - (d/2)) 



Proof: The proof proceeds analogously to the proof of Theorem 13.71 Using the full expansion, 
Lemma 13.21 instead of Lemma 12.241 leads to the following generalization of (|6.20p : 

ficontrib(w) = Z-'l^) ^ /^e-"-Mm,^,n)y-g(y)-^-^n/i,(y)-"^-"dy (6.23) 
+0{\z-''\\r\-''+^'-'') 

where the sum is over the finitely many terms with \m\ — h£ — kn < N and terms with \m.\~h£ — kn > 
N' are seen by the big-0 lemma to contribute 0(|z|^'"|r|^'^+^'*^^ ). Again, changing variables to 
y = iy' shows the summand to be equal to 



Z-''il'"l-|Ple-*^"c(m,^,n)J^-i 



The first conclusion of Proposition 16 . 21 identifies this inverse Fourier transform as an iterated deriva- 
tive (introducing a factor of i''"!), hence the summand becomes 

f)m ( 1 \ 

Z-"e-""i-|Pl(-l)^+l'"lc(m,£,n)— . (6.24) 

From Corollarv l6.5l we have 

(q-^-^) = e-i-lY |M|A*(r)-+^-W^) 

j e ^ 22(s+t)-i^{d-2)/2Y{^s + l)T{s + i+l--{dl2))- ^ ' 

Multiplying the numerator and denominator of by 11^=1 and applying once more the 

first part of Proposition 16.21 we see that 
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proving the proposition. 



□ 



The rest of the work is in determining h from the derivative (|6.22p . Begin with two classical 
regularity lemmas. 

Lemma 6.7. Let f G L^{iA + Then J'^^itAf) is standard and locally Lipschitz. 
Proof: Let g have compact support in (M'')*. 

:= (2^)-'^ / /(x)g(x)dx 



{2^)-' /(x) 



9ir) 



(2^) 



. f We" 
/(x)e- 



■ dr 



dx 
dr 



by Fubini's theorem, since e"^'''^ is bounded as r varies over the support of g and the imaginary 
part of X varies over any bounded subset of A. This shows that J-"^{iAf) is the standard function 
ih where h{r) = {2tt)~'^ J^^ /(x)e^''''^ c?x. To check the local Lipschitz condition on h, note that 



|Mr)-Mr')l = (2^)-' / /(x)(e 

Jiri+R'' ^ 

If r, r' vary over a compact set K , and x. = irj + ^ then there is a bound independent of ^: 
which implies \h{r) - h{r')\ < Cr ■ ||/||i ■ |r - r'|. 



□ 



We also require the Paley- Wiener Theorem, stated as ABGTOl . Theorem 2.5]. A generalized 
function is said to have support in a closed set K if it annihilates test functions vanishing off of K. 

Lemma 6.8 (Paley- Wiener Theorem). Suppose that A contains the convex cone K. Then the 
support of J-^^{LAf) is contained in the negative dual cone —K*. □ 



With these in hand, let K* be a connected component of the non-obstructed subset of N*. Fix 
r G K*. Suppose that for each I < j < k, there is a line segment {r + XLj : A G [0, A*]} (where 
A* could be negative) such that r -I- A*Lj is on the boundary of N* and r -|- j3\*Lj is in K* for 
all < /? < 1. In other words, for each j, traveling from r in the directions ±-Zjj, we come to the 
boundary of K* at the same time as we come to the boundary of N*. Define the integral operator 
Ij on functions on r-space by 

Ii{g)[v) = j\{v + XL,)d\. 
Let IP denote the composition over j of powers I^' . We then have the following result. 
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Theorem 6.9. Under the above geometric conditions on the component K* of the non- obstructed 
set o/N*, if 2s > d+1, then h{r) in (|6.21[) is given by 

i-\p\ IP [T-\q-^-^)] 
where p = n + nl as in Provosition 1 6. 61 

Proof: Under the condition 2s > d + 1, one has [^(x)"^! = 0{\x\~'^~'^) for some e > 0. Hence 
the function (j"" is integrable away from its poles, as is therefore 11^=1 Moving to 

ir] + W^, we avoid all poles and hence q~''~^ G L^. By Lemma [6.71 the Fourier transform is a 
standard, locally Lipschitz function, hence continuous. The domain of analyticity of / contains every 
cone whose closure is in the interior of tanx„^^^(i?). By the Paley- Wiener Theorem, therefore, the 
inverse Fourier transform vanishes outside of the closed negative dual cone, N*. Each differential 
operator d/dLj in (|6.22p may now be inverted uniquely, due to the boundary condition of vanishing 
at r + X^Lj. The unique inverse is Ij. Together with (|6.22p . this proves the theorem. □. 

Remark. Without the hypothesis 2s > d + 1 we still have 

where Ij are anti-derivative operators whose boundary conditions are not determined by continuity 
from the Paley- Wiener Theorem. 



6.6 Proof of Theorem 13.91 

We have seen that a linear factor L(x) in the denominator corresponds to a convolution with a 
Heaviside function, or equivalently, to an integral operator 1^. The integrability hypothesis in 
this result is unfortunately somewhat restrictive, ruling out, for example, the case s = l,d = 3. 
Moreover, from a computational viewpoint, it is not desirable to have the answer represented as an 
(iterated) integral. It is ther efore wor th exploring a general method for reducing the dimension of 



the integral in question. In [ABG70|, homogeneity of the integrand is exploited: integrating out 



the radial part, the Fourier transform is reduced to an integral over a cycle in (d — l)-dimensional 
projective space, which will be either the Leray cycle or the Petrovsky cycle. To this device, we add 
a residue computation that further reduces the dimension by one. Evaluation of the resulting one- 
dimensional integral leads to Theorem 13. 91 Computations in projective space rely on some standard 
constructions and notational conventions which we now introduce. 

Let TT : \ {0} CP**^^ be the projection map. Any meromorphic form lo on CP^^^ pul ls back 
to a form on \ {0}. The puUback tt* is one to one onto its range. It is well known e.g.. 



GH7g 



page 409], that the range is the set of all meromorphic [d — l)-forms on whose contraction with 
the Euler vector field ^ Xid/dxi is zero and that are homogenenous of degree zero. Here, the degree 
of fdxj^ A • • • A dxj^ is deg/ -I- k and for (d — l)-forms, those forms killing the Euler field are in 
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a one-dimensional subspace of the {d — l)-dimcnsional cotangent space at each point. Forms on 
£ipd-i jj^g^yg natu ral name s of their own, so we name them by identifying with their puUbacks 
to C*, as is done in ABG70l | and elsewhere. For computational purposes, when integrating over a 
chain C in CP''"^, we usually use an elementary chart map from a slice of C*, such as tt restricted 
to (zi, . . . , Zd~i, 1). If fj are homogeneous of degree 1 — d, for example, pulling back by this chart 
map yields 



/ fj (2) dzi A ■ ■ ■ A dz-j ^ ■ ■ ■ ^ dzd = / 
Jc Jc 



fddzi A - ■■ A dzd-i 



where C is the (unique) lifting of C to the (simply connected) slice. 

The proof of Theorem 13.91 begins analogously to the proof of Theorem 13.71 Use the general 
expansion in Lemma 13.21 just for the leading term, to write 

P(xmax + iW + y) _ p{z) ^ 



max 

+ «w + y) q{y)h{y) 



where R = 0{\y\ ^) on Cw Use Lemma [6.11 to replace Cw by Cq and use the big-0 lemma on the 
integral of R to see that 

ficontrib(w) = ^ ^.[f^ f exp(-r-y) I dy +Oi\r\-'). 

Changing variables by y = iy' and noting that dy /[q{y)h{y)] = dy' /[q{y')h{y')] gives 
6contrib(vif) .s. / exp(-zr ■ y) dy+0(|rp^) 

= Z-^P{Z)i-^T-^(X\+0{\r\-^). 



.qh 

Comparing to (|3.9|) . it suffices, therefore, to show that 



J'-' = arctan . (6.26) 

\qhj TT \ h*ir) ' 



The Leray and Petrovsky cycles 

We are left to compute the inverse Fourier transform of l/{qh). The first part of this computation, 
reducing to the Leray cycle, is valid for any hyperbolic polynoimal in any non-obstructed direction, 
so we do it in this generality. Let P/H be the ratio of two homogeneous polynomials and assume H 
is hyperbolic. Later we will specialize to the case where H = qh. Denote by :— deg H — deg P — d 
the inverse degree of homogeneity of the form {P/H) dz. 
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Let B he a cone of hyperbolicity for H and fix u G B. Fix a non-obstructed vector r G N* := 
—B*. Recall from our homotopy constuctions that there is a vector field 77 on M.'^ with the following 
properties. 



1. ?7 is homogeneous of degree +1. 

2. There is a 1-homogeneous homotopy {rjt : < t < 1} between rjo = u and 771 — r/ such that for 
all t and all nonzero y, H{iy + rjt{y) 7^ 0. 

3. For aU 2/ 7^ 0, r • -qiy) = 0. 



Indeed, a similar homotopy with the third condition replaced by r • ?7(y) < is constructed in 
Section [5] Stopping the homotopy at the instant, depending on y, that it crosses the hyperplane 
orthogonal to r, yields the desired 77 along with a homotopy as prescribed. 

Let S+ denote the hemisphere {y G K'' : |y| = 1, r • y > 0}. Let 5_ denote the other hemishpere, 
where r • y < 0. Let cycles a± be (singular triangulations of) S± oriented in such a way that 
(9(cr+ + = 0. Then a := (7+ — (t_ is a (d — l)-chain supported on S"^^^ whose boundary is 
supported on the equator {y G S"^^^ : r • y = 0}. 

The map </> defined by (/)(y) := iy + viv) induces a covariant map 0, on cycles and homology. The 
chain (f)*{a) maps to a cycle in with boundary in the complex hyperplane := {z : r • z = 0}. 
Hence (/)*((t) represents a homology class in iJd_i(C'', Xr). For any homogeneous set W C C*, let 
W denote the projection ttW of W to CP''"^. The sets Vq and are homogenous, therefore the 
pair (C'' - Vq, (C^ - g) n X^) projects radially under tt to {Cf^-^ - V^, (CP'^-i - W) ^ ^)E1 

Definition 6.10 (Leray and Petrovsky cycles). The chain a — oi{r) :~ TT(j){a) is called the Leray 
cycle and its class tt^^^ct] G H'^-'^{CF'^-^ - V^, (CP'^-i - V^) n ^7 ) is called the Leray class. 

The boundary of a is a cycle (3 representing a class in Hd-2{'CP'^'^^ — Vq) n X^). Define the 
Petrovsky cycle ^ to be a tubular neighborhood around (3 orthogonal to Xr- This is the image under 
of a cycle supported on {y : |r • y| = e} and avoiding Vq. 



The following result is proved in ABGTCt Theorem 7.16]; here we correct a typo: the second 



ap pearance of namely the one in (7.17'), should be just Xg! see equation (1.6) on page 122 



of 



d, by 



ABG70|. Define a (d — l)-form killing the Euler vector field and having homogeneous degree 



- ^(-l)-'+^Zj dzi A • • • A dzj A • • • A dzd . 

.7=1 



*The bars denoting projective varities are about to proliferate; we apologize for the mess, but we tried dropping 
them but we became confused about which varieties were projective and which were afBne. 
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To explain what is about to appear in (|6.27p - (|6.28p below, we must see why the integral of a 
meromorphic 2-form on CP^ over a relative homology class in i72(CP^, A^) is well defined. Indeed, 
integration over a relative homology class with respect of a complex submanifold of positive co- 
dimension is always well defined, for the following reason. Let C be a representating cycle for the 
class, that is a 2-chain with boundary in A4. The integral over any other relative cycle differs from 
the integral over C by the integral over a relative boundary, a relative boundary being an absolute 
boundary plus something in Ai. Since duj — for any meromorphic 2-form on CP^, the integral 
over the boundary vanishes, and since M has positive complex co-dimension, the second part of the 
integral vanishes as well. 

Theorem 6.11 (Reducing Fourier integrals to Leray/Petrovsky cycles). Let P/H he hyperbolic and 
fix u in a cone B of hyperbolicity for H as above. Let be the inverse degree of homogeneity of 
the form {P/ H)uj, that is, — deg H — deg P ~ d. Let a be the Leray cycle and 7 be the Petrovsky 
cycle. If d^ > then 

/ P\ jd, + l I- p 

(77) -(2iF^ 77- ^'-''^ 

□ 

Remarks, (i) Note that the introduction of the factor (r • z^' makes the integrand 0-homogeneous, 
which is exactly what we need to interpret it as a form on CP''^^. (ii) In the case d* < 0, the 
integrand contains a negative power of r • z. Let Res(ci;) be |d*|*'* residue of the integrand (r • 
•z.Y'{P/H)lo along the projective hyperplane Xj.. The product structure in the Petrovsky cycle 
immediately reduces the integral one dimension further to Res((jj). In the case d* > 0, the integral 
does not localize to the boundary cycle (3 and one must work harder to kill one more dimension. 



while if d^ < then 



Residue reduction 

The second step is to reduce by one further dimension via a residue computation. The first half of 
this step still works in any dimension. Begin by observing: 

Lemma 6.12. The homology group Hd-iiCF'^-^Xr) vanishes. 

Proof: CP'^^-^ n Xj. is homeomorphic to CP*^^^. If p < 2(d — 2) then this inclusion induces an 
isomorphism Hp{CP'^^^ n X^) Hp{CP'^^^), where both groups have rank 1 if p is even and vanish 
otherwise. It follows that the first and last arrows are isomorphisms in the exact sequence 

ifd-i(CP'^~' nXr) ^ Hd-iiCF''-') ^ Hd-iiCF''-\Xr) ^ Hd-2{CF''-^ n Xr) ^ i/d-2(CP'^"l) 

hence the middle group vanishes. □ 
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We now make use of the Thorn isomorphism to "pass a through Vh" and obtain a {d~ 2)-chain, 
whose tubular neighborhood is homologous to a. The details are as follows. By Lemma [6.121 any 
chain representing the pLeray class is a boundary of a d-chain (3 G {CW"^~^ , X^)- We may therefore 
choose a d-chain C in CP''^^ whose boundary is a plus something in Xj- (because dC is a cycle and 
a is not, the part in X^ will be nonzero). Perturbing C if necessary, we can assume that C intersects 
Vh transversely. The dimension of the intersection of the c?-chain C" with the surface Vh having 
co-dimension 2 is a (d — 2)-chain 6, whose orientation is prescribed by the orientations of C, CP^~^ 
and Vh- The chain S has boundary in Vh H Xr, and is therefore a relative cycle in {Vh, Vh H X^)- 

Now we are at a point where we require the dimension to be 3. The chain C has (real) dimension 3 
and the surface Vh has dimension 4. The projective variety Vh may not be smooth, but its singular 
set has complex co-dimension at least 1, hence has real dimension at most 2. Generically perturbed 
C therefore does not intersect the singular set of Vh, and hence S is supported on the set of smooth 
points of Vh- We may define the tubular neighborhood T{S) supported on the set {\H\ = e}, which 
is locally a product of S with a small circle about the origin in with the standard orientation. 

Integration around this circle is computed by taking a residue. We recall a definition of the 
residue form on any complex space. Let 9 be a meromorphic form with a pole on the set Vh, the 
pole being simple except on a proper subvariety Vh H Vk- Then the residue is defined as follows. 
Write 9 = H ■ uj where uj is holomorphic away from on Vk- Define Res[0, V^] to be the unique form 
on Vh that satisfies 

Res[9,\^]AdH = uj 

away from Vk- In coordinates, the residue of {G/H)dzi A • • • A dzd is given by [G / [dH / dzi))dz2 A 
• • • A dzd; (with Zj in place of zi, we have the alternative expression {~iy^^{G/dH/dzj)dzi A • • • A 
dzj^i A dzj^i A • • • A dzd)- The following well known result may be demonstrated by expressing the 
integral over the tube as an iterated integral, first around a circle. 

Lemma 6.13. Suppose d = 3. The Leray class a is homologous to the tube T(S) around S. Conse- 
quently, for any meromorphic form 9 on H with a simple pole at H , 



9 = {2ni) [ Res[9;VH] - 
Js 



a •'T(S) Js 

In particular, when d* = m Theorem \6-lll putting this together with (|6.27p and specializing to 
H = qh and 9 ~ lo / {qh) yields 



where 



1 \ i f UJ — 1 



qh) (27r)2 qh 27r 

■= Res J; 



□ 
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Remark. In higher dimensions, it may happen that S intersects the singular set of V^- In that case, 
one might expect a version of Lemma |6 . 1 31 showing the Leray cycle to be homologous to the sum of 
a tubular neighborhood of d away from the singular set and a cycle supported on a neighborhood of 
the singular set. 

The special case of a cone and a plane 

For the remainder of this section, we specialize to the case in Theorem 13.91 In addition to d = 3, 
we suppose k = s = 1, whence the inverse degree, c?», of homogeneity is zero. We further specialize 
to the case where Vg fl D R"^ ^ and where the variety V-.^^ has precisely two points, each of 
multiplicity one. Thus we are in the case of Section [5.31 which arises in several of our applications 
and is illustrated in figure [141 The normal cone is shaped as in figure 13.31 and the vector field rji 
constructed in Section [5.31 is one we can use to construct the Leray cycle. In this case {d — 3), since 
d avoids V^, we may restrict to Vq and write 



Let p be a common zero of q and h. The space of 2-forms on CP^ is one-dimensional, hence at 
p, the form w is a multiple of the form dq /\ dh. The double residue 



where and are the values of the parameter t for which h vanishes. 

Proof: The form is meromorphic on Vq with precisely two simple poles. Therefore, it may be 
written as C[dt/{t — ts) — dt/{t — ^4)]. Taking the residue at t = ^3 yields Res{uJL',t3) = C. But 



COL := Res — ; V, 



Res'2' :=Res — ,V,~n 

.qh 

is the value of this ratio, that is, uj = Res^^-* (p) dqhdh at p. We have the following explicit description 
of the residue form lj^. 



Lemma 6.14. Let t : CP^ 



Vq be any local parametrization. Then 
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iterated resides are the same as multiple residues, hence 



= Res 

= Res^^ 
= Res(2). 



qn 



To carry out the rest of the computation, let us choose coordinates for in which q 



-y 



Although these are unrelated to the coordinates in which h and r are described, we will use them 
to compute a coordinate-free description of the integral. 

The cones of hyperbolicity of q are the two components of > + z^, one containing the 
negative x-axis and one containing the positive x-axis. Recall that the cones of hyperbolicity for qh 
are these cones, bisected by the plane h = 0. Recall we have fixed u e i?, where B is one of these 
sliced cones and its dual has a teardrop shape. 



The space Vq is a quadratic curve in ' 



We may choose the explicit parametrization t 



Vg by 

{x,tx) ^ {f^ + l,2t,t2 _ 1) 

in the (a;, y, z)-coordinate system. Topologically, CP^ is a sphere with WS'^ as its equator, and our 
parametrization has the nice feature that the copy of RP^ inside CP^ maps to the real part of the 
quadric Vq. Thus figure [19] depicts Vq as a sphere and the real part as the equator. The points of Vq 
where r • x vanishes arc denote ti and t2 and the points where h vanishes are denoted by and ^4. 
The assumption that Vj^ and Vq intersect in real space is equivalent to h lying outside the dual cone, 
which is equivalent to t having real roots. Thus ^3 and are shown on the equator in figure 1191 
Note that the two arcs into which ts, ^4 separate the equator differ in their position visavis the cone 
of hyperbolicity: one of them bounds it (we will call this arc the active one), while the other not. 

If r is in the dual cone then ti and ^2 will be complex conjugates, while if r is in the pointy region 
of the teardrop then ti and ^2 will be real. 



Case 1: ti and t2 are complex 

We are nearly ready to evaluate the integral, but we need first to understand the cycle 5. The 
intersection class ^ is a relative cycle in (Vg, Vq n X^)- Thus we may draw a representative of 
this class as a path, beginning and ending in the set {ti,t2}- The meromorphic residue form lol 
is holomorphic away from and ^4 where it has simple poles. The integral of this form over 5 is 
therefore determined by combinatorial invariants of 5: the positions of the endpoints and the number 
of signed intersections with the two equatorial arcs bounded by and t^. 
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Figure 19: the topological sphere Vq, its real part (the equator) and its intersections with the planes 
{h = 0} and Ji^ 

Lemma 6.15. The homology class of S in {Vq \ {ta, t^}^ {ti, 12}) is that of an oriented path from t2 
to ti, intersecting one equatorial arc exactly once. 



Proof: It is shown in ABG7Cll | (see in particular figure 6b there and the paragraph preceding it) 
that one can find a representative of the Leray class such that its boundary (in — Vq O Xr) is 
localized near the complex points ofVq D X^, i.e. in our situation is the sum of small circles around 
the complex zeros of q in Xr, oriented according to their imaginary parts (note that h has no non-real 
zeros there). It follows that the boundary of the (relative) cycle 6 is given by 

dS = [ti] - [t2] 

where ti has the positive imaginary part, and 5 is the claimed path, plus one or more absolute cycles 
(i.e. oriented closed loops) in Vq and V^. 

To find the homology classes represented by these arcs and loops, we recall the definition of 
the Leray class: the vector field constructed at the end of section 15.31 is restricted to the unit 
sphere, defining the relative class, the sum of oppositely oriented semispheres separated by the 
hyperplane Xj.. To evaluate 5 we find a homotopy shrinking these hemispheres to a point, keeping 
their boundaries in Xj. and tracking where the resulting 3-chain hits Vq and V^. This homotopy 
proceeds in two stages: first we take the linear homotopy of (the restriction to the unit sphere of) 
77, the vector field constructed in Section [?751 to the (restriction to the unit sphere of the) constant 
vector field x defined in the same place. In the second stage we collapse the sphere to a point, 
keeping the constant vector field. 

We note first that in this deformation yields the empty set: at no instant arc the deformed 
vectors tangent to h, which would be necessary if the deformation were to intersect h. This is not 
the case for q, and indeed, we know already that the boundary of 5 there is nontrivial. The class of 
5 is completely determined by the index of intersection with the active arc between t^ and ^4. The 
intersection number of 5 with the active arc is just the number of points in the real part of Vq where 
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point of vanishing of tlie vector field 



Figure 20: showing where the homotopy intersects the quadric 



our deformation results, at some instant of the homotopy, in a vector field tangent to the (real part 
of) the quadric. It is immediate that there is a single point and a single time in the homotopy where 
this occurs (in fact, the vector field vanishes at this place and time), and it is easy to check that this 
uniquness survives small perturbations (see figure [^D|) . □ 

Having this geometric understanding of Vq, ti, t2, ^3, ^4 and d, we may now compute the integral. 
We find that 



/ Res(2) . ( ^ Res^^J loe 

JS V-h t-Uj ■ ^ (il -i3)(<2 -t4) 



-,(2) - ^3)(^1 - tj) 



^ i{a + P) Res^^^ . 

Here, the fact that ii, t2 are complex conjugates while <4 are real implies that the numerator and 




Figure 21: the logarithm as an arctangent 



denominator are complex conjugates and the logarithm is purely imaginary, being in fact 2i times 
the arctangent of the sum a + f3 oi the angles shown in figure [2TJ The logarithm is therefore given 
by twice the arctangent of the ratio of imaginary to real parts in the numerator. 

Whenever ti,t2 satisfy a quadratic equation + at + b with real coefficients while ^3, ^4 satisfy a 
quadratic equation t'^ + a't + 6' = 0, also with real coefficients, then simple algebra shows the cross 
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ratio to be given by 



ih - h)(t2 ~ ts) ^ b + b' - aa'/2 + l^/^^^/Ah' - a'^ 
{ti-h){t2-U) b + b' - ad 12 - i^a^ - 4&V4fo' - a'2 ' 

The ratio of the imaginary to real parts simplifies considerably, so we obtain the equivalent expres- 
sions 

2i arctan — % ^ — , . (6.30) 

6 + 5' - aa'/2 ^ ' 



Here, we recall the definition of the range of the arctangent function in Theorem 13. 9i namely < 
arctan a; < tt. 

Let L = £ix + i2y + (sz describe h in our coordinate system. Then ia, solve q — L = 0. The 
minimal polynomial for ^3 and t4 (produced, for example, in Maple as an elimination polynomial for 
the ideal {x — (t'^ + l),y — 2t, z — {t^ — 1) , x"^ — y"^ — , lix + €2!/, f-?,z) is given by 

f H —t + — i . 

Similarly, let r = r\x + r^y + ^32, giving the minimal polynomial for t\ and as 

9 2r2 r\ — r3 
+ —t + — . 

n +r3 ri+ ra 

Plugging in a = 2r2/(r3 + n), 6 = (n - r3)/{n + rs), a' = 2^2/(4 + d) and b' = - £3)7(^1 + £3) 
to (|6.30p now gives 



LOL = 2iRes(2) arctan 2 3V 12^3 

y rill - r2i2 - rzl^ j 

where the two quantities under the radical signs are both positive. Writing the right-hand side as a 
combination of coordinate-free quantities, this becomes 



V'g*(r,r)J-g*(/i, h) 

LUL = 2i Res^^' arctan ^— = . (6.31) 

g*(r,/i) 



Combining this with the result of Lemma 16.131 shows that 



T-^ [ ^\ ^ -Res(^) arctan .V^V-^-(^^ 



,q-hj TT q*{r,h) 
and checking this against [6.261 proves Theorem 13.91 in the case where r is inside the dual cone. 



Case 2: ti and t2 are real 



In this case, r is in the pointy region of the teardrop. The result (6.23) from ABG7Cll | tells us 
that the Leray cycle, which is by definition a relative cycle, is an absolute cycle. It follows that 
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the intersection class 5 is an absolute cycle in the twice punctured sphere. Thus ti = t2 when the 
roots are real, and the cycle S is represented by a circle. Geometrically, in the previous case (ii,i2 
complex), as r approaches the boundary of the dual to the conic, the points ti and ^2 converge to a 
single point on the equator and the arc S closes up into a circle. Provided this point of convergence 
is not one of the poles, ^3 or t^, the integral will approach a well defined limit, which is the integral 
over the absolute intersection cycle. By continuity, the homology class of this absolute cycle cannot 
vary as the limit point varies over the common boundary of the two regions of the teardrop, nor can 
this class vary as r varies over the pointy region of the teardrop. 

To summarize, there is a constant c such that for all r in the pointy region, the integral is 
equal to c. This is also the limiting value if r approaches any point of the common boundary from 
inside the other region, and thus concides with the limit of the quantity in the previous case, as r 
approaches any ray in the common boundary; in the limit the arctangent is tt and we obtain simply 
P(Z)Res(2^Z-"-. □ 

Remark. If r crosses out of the dual conic at a point a not on the boudary of the pointy region, 
it exits the normal cone. We know the integral must become zero in this case. Geometrically, this 
corresponds to ti and ^2 coming together in a cycle homologous to zero. There is a discontinuity 
if a is one of the two projective points of tangency in figure [101 Near these two points, ti and ^2 
approach ^3 and respectively. Crossing out out of the dual cone on one side or the other will 
cause S to close up to a null or non-null cycle, in the former case the integral is zero; the difference 
between the two integrals is the residue at the pole or t^. 



7 Further questions 

1. It would be nice to remove the integrability hypothesis 2s > d+1 from Theorem 16. 91 Doing so 
would necessitate a specification of which antiderivative 5(r + Aij) dX is meant when the 
integral in question is not convergent. There are cases when the is one "obvious" interpretation 
of this as a closed form function, but proving this to be correct requires better understanding 
of the generalized function partially identified as g on the Paley- Wiener cone. 

2. When d — 3 but Q has an isolated singularity of degree greater than 2, the techniques of 
Section [6.61 are still applicable through Lemma l6.13l The representation in Lemma 16.141 must 
be replaced by one with four poles, and the intersection class in Lemma |6. 151 correspondingly 
specified. Work is in progress on details of this computation and its applicatoin to the Fortress 
tiling ensemble. 

3. In principle, generalized Fourier transform theory should give us some information on asymp- 
totics in scaling windows near the boundary or in obscured directions. An additional complica- 
tion is that, because projective homotopies do not exist giving exponential decay of exp(— r-x) 
for these r, one must verify the existence of chains on which exp(— r • x) decays sufficiently 
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rapidly to justify the exchanges of hmits in Lemma 16.31 and elsewhere. This, together with 
the increased complexity of Fourier integrals with varying parameters, has kept us thus far 
from obtaining limit theorems near these boundaries. This is perhaps the most broad and 
challenging open problem pertaining to the results in this paper. 
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